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INTRODUCTION 


To  bring  our  reporting  cycle  into  step  with  the  new  JSEP  schedule,  this 
report  covers  progress  made  in  only  the  9-month  period  1  April,  1985  to  31 
December,  1985  for  the  eleven  research  units  funded  under  the  Joint  Services  Elec¬ 
tronics  Program  at  Harvard  University.  It  is  broken  down  into  four  major  divi¬ 
sions  of  electronic  research —solid  state  electronics,  quantum  electronics,  informa¬ 
tion  electronics,  and  electromagnetic  phenomena.  It  also  includes  a  report  of 
Significant  Accomplishments  which  contains  selected  noteworthy  results.  These 
are:  “Ultra-High-Speed  Real-Time  Optical  Processing  Via  Convolution  and 

Cross-Correlation  in  Inhomogeneously  Broadened  Optical  Materials”  by  Research 
Unit  6  and  “Electromagnetic  Missiles”  by  Research  Unit  11. 
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I.  SOLID  STATE  ELECTRONICS 
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1.1  Electronic  Structure  and  Magnetic  Properties  of  Diluted  Magnetic 
Semiconductors  (DMS).  H.  Ehrenreich,  K.C.  Hass,  B.  Larson,  and  R. 
Lempert,  Contract  N00014-84-K-0465;  Research  Unit  1. 


The  group  has  been  investigating  hypothetical  zinc-blende  MnTe  in  the  Type 
I  and  Type  II  anti-ferromagnetic  states  as  well  as  the  disordered  magnetic  state. 
These  calculations  are  complemented  by  ASW  first  principles  computations. 

These  results  serve  as  an  input  for  careful  empirical  tight-binding  calculations  for 
disordered  CdMnTe  and  permit  investigation  of  the  pressure-dependence  of  the 
electronic  levels  and  the  cohesive  energy  of  the  alloy.  They  also  serve  as  a  basis 
tor  the  investigation  of  the  magnetic  properties  as  described  below. 


The  fundamental  exchange  constants  responsible  for  the  magnetic  properties 
of  diluted  magnetic  semiconductors  are  also  being  investigated  [l] .  The  d-d 
exchange  interaction  has  been  shown  to  be  dominated  by  superexchange  and  to  be 
far  more  important  than  other  exchange  interactions.  The  strength  of  the  sp-d 
exchange  interaction,  which  determines  the  optical  properties  of  practical  interest 
in  the  band-gap  region  and  the  properties  of  magnetic  polarons  has  also  been 
determined.  The  results  appear  to  be  in  good  agreement  with  experiment. 

Reference: 

1.  B.E.  Larson,  K.C.  Hass,  H.  Ehrenreich,  and  A.E.  Carlsson,  “Exchange 

Mechanisms  in  Diluted  Magnetic  Semiconductors,”  Solid  State  Commun  56. 
4,  347  (1985). 

1.2  Magnetization  Steps  in  CdMnTe.  K.C.  Hass,  B.E.  Larson,  and  H. 
Ehrenreich,  Contract  N00014-84-K-0465;  Research  Unit  1. 

More  recently,  the  magnetization  steps  in  CdMnTe  observed  at  the  National 
Magnet  Lab.  have  been  interpreted  [l ] .  These  are  exhibited  in  measurements  of' 
the  magnetization  as  a  function  of  magnetic  field  at  very  high  field  of  R.L. 
Aggarwal  and  result  from  spin  flips  modified  by  nearest-neighbor  exchange 
interactions.  Magnetization  steps  have  been  observed  in  alloys  having  low  Mn 
concentrations.  Our  calculations  are  being  extended  to  higher  concentrations 
using  CPA  techniques.  The  experimental  observations  require  higher  fields. 
Despite  experimental  difficulties,  however,  it  is  important  to  understand  these 
effects  in  order  to  make  a  quantitative  determination  of  the  strength  of  the  inter¬ 
nal  magnetic  fields  which  affect  both  transport  and  optical  properties. 

Reference: 

1.  B.E.  Larson,  K.C.  Hass,  and  R.L.  Aggarwal,  “Effects  of  Internal  Exchange 
Fields  on  Magnetization  Steps  in  Diluted  Magnetic  Semiconductors," 
(accepted  for  publication  in  Phy.s  Rev.  B). 
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1.3  Semiconducting  Pseudobinary  Alloys.  H.  Ehrenreich,  R.  Lempert,  and 
K.C.  Hass,  Contract  N00014-84-K-0465;  Research  Unit  1. 

The  MCPA  has  been  applied  to  the  anionic  II- VI  alloy  Zn(Se,Te)  and  the 
cationic  III-V  alloy  (In,  Ga)As  [l].  Detailed  comparisons  have  been  made  of  the 
energy  levels  of  chalcopyrite  and  corresponding  disordered  structure  for  the  II- VI 
compound  in  order  to  investigate  the  possible  effects  of  short-range  order.  The 
ordered  crystal  approximation  is  not  expected  to  be  valid  in  the  understanding  of 
band  edge  properties  since  effects  such  as  disorder-induced  bowing  and  disorder 
scattering  on  the  electron  mobility  are  neglected.  These  investigations  have  also 
shown  that  the  MCPA  is  a  better  approximation  to  the  alloy  band  structure  than 
the  bond-centered  CPA  used  by  Chen  and  Sher. 

As  pointed  out  in  last  year’s  comprehensive  proposal,  the  MCPA  is  suitable 
for  treating  both  atomic  and  chemical  disorder  in  pseudobinary  alloys  such  as 
Ga(In,  As).  Its  use  was  motivated  in  the  present  case  by  the  experimental  obser¬ 
vation  that  in  this  alloy  the  GaAs  distances  and  the  InAs  distances  correspond 
closely  to  those  encountered  in  the  pure  crystal,  b”t  that  the  average  lattice  con¬ 
stant,  as  determined  for  example  by  x-ray  measurements,  closely  follows  Vegard’s 
Law.  In  addition  to  the  combined  effects  of  chemical  and  structural  disorder  on 
electronic  structure,  Ehrenreich  and  Hirth  [2]  have  also  explained  the  experimen¬ 
tally  observed  dislocation  reduction  in  GaAs  upon  the  addition  of  small  amounts 
of  In  using  a  classical  strain  model.  The  work  pertaining  to  charge  dislocations 
and  jogs  in  HgCdTe  described  in  last  year’s  proposal  has  been  completed  and  pub¬ 
lished  [3]. 

References: 

1.  R.  Lempert,  “Electric  Theory  of  Semiconducting  III-V  and  II- VI  Alloys,” 
Ph.D.  Thesis,  Harvard  University,  Division  of  Applied  Sciences,  1985. 

2.  H.  Ehrenreich  and  J.P.  Hirth,  “Mechanism  for  Dislocation  Density  Reduction 
in  GaAs  Crystals  by  Indium  Addition,”  Appl.  Phys.  Lett.  46,  7,  668  (1985). 
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3.  J.P.  Hirth  and  H.  Ehrenreich,  “Charged  Dislocations  and  Jogs  in  Hg^Cd^Te 
and  Other  II- VI  Compounds,”  J.  Vac.  Sci.  Technol.  A3,  2,  367  (1985). 

1.4  Electronic  Structure  of  HgCdTe  and  Related  Non-Magnetic  II- VI 
Alloys:  HgZnTe,  CdZnTe.  K.C.  Hass,  Contract  N00014-84-K-0465; 
Research  Unit  1. 

K.C.  Hass  is  investigating  the  stability  of  ordered  50-50  alloys  of  these 
materials  having  the  chalcopyrite  structure  and  the  corresponding  disordered  case. 
This  investigation  is  important  because  the  influence  of  Cd  d-levels  on  the 
cohesive  energy  is  significant  and  is  not  yet  understood.  Furthermore,  it  permits 
an  investigation  of  the  possible  importance  of  short-range  order  among  the  consti¬ 
tuents  as  contrasted  to  possible  segregation.  It  may  be,  for  example,  that  the  50- 
50  alloy  is  less  stable  than  the  phase  separated  mixture.  Charge  transfer  effects 
which  are  important  because  they  stabilize  the  cohesive  energy  can  be  investi¬ 
gated.  These  transfers  have  marked  effects  on  electronic  levels  and  influence  band 
offsets,  which  are  of  particular  importance  in  heterojunctions.  Finally,  these  cal¬ 
culations  serve  as  another  input  to  the  parameterized  alloy  calculations  described 
above. 

1.5  Chemical  Trends  in  Diluted  Magnetic  Semiconductors.  H.  Ehren¬ 
reich,  B.  Larson,  and  K.C.  Hass,  Contract  N00014-84-K-0465;  Research  Unit 
I. 

We  are  also  investigating  the  chemical  trends  in  the  family  of  materials 
characterized  by  CdMnTe  and  ZnMnTe  with  Se  and  S  substituted  for  the  chal- 
cogenide.  We  have  been  able  to  explore  the  variation  of  the  position  of  the  d- 
levels  and  the  strength  of  the  sp-d  hybridization  by  utilizing  simple  but  rather 
generally  applicable  results  obtained  from  detailed  calculations  for  CdMnTe.  The 
same  applies  to  the  magnitude  of  exchange  interactions,  which  are  also  under 
investigation. 


1.6  Amorphous  Semiconductors.  H.  Ehrenreich  and  K.C.  Hass,  Contract 
N00014-84-K-0465;  Research  Unit  1. 

Sum  rules  involving  the  dielectric  function  weighted  by  various  powers  of  the 
frequency  have  been  utilized  to  analyze  the  differences  between  chemical  bonding 
in  amorphous  materials  and  their  crystalline  counterparts  [l] .  They  are  seen  to 
provide  fairly  direct  information  concerning  the  short-range  order  present  in 
amorphous  semiconductors.  They  should  therefore  be  useful  as  a  method  of 
experimental  investigation  of  such  effects.  The  so-called  independent  band  model 
for  amorphous  semiconductors,  which  supposes  that  the  valence  and  conduction 
bands  are  independent  of  each  other  is  frequently  invoked  in  the  interpretation  of 
the  optical  fundamental  absorption  edge.  Tauc’s  law  is  based  on,  and  Urbach 
Tails  have  been  explained  using  this  model. 

Reference: 

1.  K.C.  Hass  and  H.  Ehrenreich,  “Electronic  Structure  Models,  Bonding  and 
Optical  Moments  in  Amorphous  and  Crystalline  Semiconductors,”  Annals  of 
Physics  164,  77  (1985). 

1.7  Photovoltaics.  H.  Ehrenreich,  Contract  N00014-84-K-0465;  Research  Unit 

1. 

Ehrenreich’s  interest  in  an  involvement  with  photovoltaics  continues,  albeit 
on  a  somewhat  reduced  level.  Last  March  he  was  asked  to  testify  before  the 
House  Committee  on  Science  and  Technology  chaired  by  Congressman  Don  Fuqua 
concerning  the  FY’86  budget  for  photovoltaics  [l]. 

Reference: 

1.  H.  Ehrenreich,  “Hearings  before  the  Subcommittee  on  Energy  Developments 
and  Applications  of  the  Committee  on  Science  and  Technology,”  U.S.  House 
of  Representatives,  March  5,  1985,  No.  31,  Vol.  III. 
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a.  Papers  Submitted  to  Refereed  Journals  (and  not  yet  published) 

K.C.  Hass,  B.E.  Larson,  H.  Ehrenreieh,  and  A.E.  Carlsson,  “Magnetic 
Interactions  in  Diluted  Magnetic  Semiconductors,”  accepted  for  publication 
in  J.  Magnetism  and  Magnetic  Materials  (Proceeding  of  the  Int.  Conf.  on 
Magnetism,  1985). 

B.E.  Larson,  K.C.  Hass,  and  R.L.  Aggarwal,  “Effects  of  Internal  Exchange 
Fields  on  Magnetization  Steps  in  Diluted  Magnetic  Semiconductors,” 
accepted  for  publication  in  Phys.  Rev.  B. 


b.  Papers  Published  in  Refereed  Journals 

H.  Ehrenreieh  and  J.P.  Hirth,  “Mechanism  for  Dislocation  Density  Reduction 
in  GaAs  Crystals  by  Iridium  Addition,”  Appl.  Phys.  Lett.  46,  7,  668  (1985). 

B.E.  Larson,  K.C.  Hass,  H.  Ehrenreieh,  and  A.E.  Carlsson,  “Exchange 
Mechanisms  in  Diluted  Magnetic  Semiconductors,”  Solid  State  Commun.  56. 
4,  347  (1985). 

J. P.  Hirth  and  H.  Ehrenreieh,  “Charged  Dislocations  and  Jogs  in  Hg[_ICd,Te 
and  Other  II- VI  Compounds,”  J.  Vac.  Sci.  Technol.  A3,  2,  367  (1985). 

K. C.  Hass  and  H.  Ehrenreieh,  “Electronic  Structure  Models,  Bonding  and 
Optical  Moments  in  Amorphous  and  Crystalline  Semiconductors,”  Annals  of 
Physics  164,  77  (1985). 

H.  Ehrenreieh,  “Hearings  before  the  Subcommittee  on  Energy  Developments 
and  Applications  of  Committee  on  Science  and  Technology,”  U.S.  House  of 
Representatives,  March  5,  1985,  No.  31,  Vol.  III. 


c.  Books  (and  sections  thereof)  Submitted  for  Publication 

H.  Ehrenreieh  and  D.  Turnbull,  eds.,  Solid  State  Physics ,  Volume  39. 
Academic  Press. 
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Books  (and  sections  thereof)  Published 

H.  Ehrenreich  and  D.  Turnbull,  eds.,  Solid  State  Physics,  Volume  38, 
Academic  Press,  1984. 


Invited  Presentations  at  Topical  or  Scientific/Technical  Society 
Conferences 

B.E.  Larson,  K.C.  Hass,  and  R.L.  Aggarwal,  “Magnetic  Interactions  in 
Diluted  Magnetic  Semiconductors,”  International  Conference  on  Magnetism 
(San  Francisco,  August  1985). 

K.C.  Hass  and  H.  Ehrenreich,  “Implications  of  Optical  Moment  Comparisons 
of  Amorphous  and  Crystalline  Semiconductors”,  American  Physical  Society 
(Baltimore,  March  1985). 

K.C.  Hass,  B.E.  Larson,  H.  Ehrenreich,  and  A.E.  Carlsson,  “Electronic  Struc¬ 
ture  and  sp-d  Exchange  in  Semimagnetic  Semiconductors”,  American  Physi¬ 
cal  Society  (Baltimore,  March  1985). 


Honors/Awards/Prizes 

H.  Ehrenreich,  Visiting  Exchange  Scholars  Program  with  People’s  Republic  of 
China,  National  Academy  of  Science,  May  i985. 

R.  Lempert,  “Electronic  Theory  of  Semiconducting  III-V  and  II- VI  Alloys,” 
Ph.D.  Thesis,  Harvard  University,  Division  of  Applied  Sciences,  1985. 

K.C.  Hass,  “Topics  in  Electronic  Theory  of  Disordered  Semiconductors”, 
Ph.D.  Thesis,  Harvard  University,  Division  of  Applied  Sciences,  October 
1984. 
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1.8  Investigation  of  Amorphous  Hydrogenated  and  Fluorinated  Alloys 
of  Silicon  and  Germanium.  J.  Hanna,  K.D.  Mackenzie,  Z.L.  Sun,  J.H. 
Burnett,  J.R.  Eggert,  Y.M.  Li,  and  W.  Paul,  Contracts  N00014-84-K-0465 
and  SERI  Subcontract  XB-2-02144-1  of  DOE  Prime  Contract  DE-AC-02-83- 
CH10093;  Research  Unit  2. 

In  the  previous  report  (No.  98),  we  discussed  our  conclusions  regarding  the 
properties  of  amorphous  hydrogenated  alloys  of  silicon  and  germanium 
(a— Si,_.,Ge,:  H).  These  alloys  have  a  smaller  band-gap  than  a-Si:  H,  which  makes 
them  useful  for  electrophotographic  materials  with  improved  “red”  response  and 
also  for  the  lower-band-gap  component  of  tandem  solar  cells.  From  the  funda¬ 
mental  point  of  view,  the  amorphous  group  4-group  4  alloys  may  be  looked  upon 
as  prototypes  for  many  compound  amorphous  semiconductors  of  the  silicon  fam¬ 
ily.  A  characteristic  of  these  4-4  (or  3-5  etc.)  alloys  is  that  their  photoelectronic 
properties  are  inferior  to  those  of  a-Si:  H,  for  reasons  that  are  not  yet  understood. 
We  had  earlier  investigated  a-Si,_,Ge,:  H  produced  by  diode  sputtering  and  d.c. 
glow  discharge,  and  this  was  followed  up  with  a  detailed  study  of  alloys  prepared 
by  the  more  conventional  technique  of  r.f.  glow  discharge.  Based  on  a  systematic 
examination  of  the  structural,  optical  and  transport  properties  of  samples 
prepared  under  a  variety  of  deposition  conditions,  we  have  published  our  conclu¬ 
sions  regarding  their  atomic  structure  and  electronic  band  structure  [l].  One  of 
these  conclusions  concerned  an  inferred  looser,  more  defected  atomic  structure  in 
the  Si-Ge  alloys  than  in  Si  itself. 

Several  years  ago,  the  Principal  Investigator  suggested  [2]  that  improvements 
in  the  properties  of  Si-Ge  alloys  might  eventuate  if  a  different  bond  terminator 
than  H  were  used.  The  bond  terminator,  F,  had  already  been  proposed  and 
investigated  by  the  group  at  LCD  [3].  This  suggestion  was  taken  up  by  a 
Japanese  group  [4]  which  determined  deposition  conditions  for  a-Si^Ge,:  H:  F 
leading  to  much  improved  photoconductive  response.  We  have  now  completed  the 
first  phase  of  an  investigation  of  amorphous  alloys  of  Si,_,Ge,  by  the  glow 


discharge  decomposition  of  mixtures  of  SiF4,  GeF4  and  H2  and  have  presented 
our  findings  [5].  In  brief,  they  confirm  the  result  of  the  Japanese  group  that  the 
photoconductivity  may  be  improved  by  about  an  order-of-magnitude,  while  they 
show  that  in  every  other  respect  (conductivity,  optical  absorption  edge  spectrum, 
vibrational  absorption  spectrum,  photoluminescence)  there  is  no  essential 
difference  between  a-Si^Ge,:  H  prepared  from  hydrides  and  a-Si^Ge,:  H:  F 
prepared  from  fluorine-containing  gases. 

The  investigation  of  a-Si^Ge,:  H:  F  from  mixtures  of  SiF4,  GeF4  and  H2  was 
carried  out  on  samples  prepared  at  different  substrate  temperatures,  with  different 
SiF4:  GeF4  gas  pressures  and  flow  rates,  and  at  different  r.f.  power  levels.  Optimum 
parameters  were  chosen,  with  a  focus  on  the  properties  of  samples  with  an  energy 
gap  in  the  1.4-1.5  eV  range.  A  wide  variety  of  properties  were  then  measured 
(vibrational  absorption,  .optical  absorption  edge,  conductivity  versus  temperature, 
photoconductivity  magnitude,  photoconductivity  spectra,  photoluminescence). 
Although  minor  differences  were  found  in  many  of  these  properties  from  those 
measured  on  alloys  prepared  from  hydrides,  the  most  striking  and  consistent 
change  was  an  order-of-magnitude  increase  in  the  photoconductivity. 

Several  points,  however,  were  notable.  First,  the  content  of  F  was  of  the 
order  of  1  atomic  percent,  while  that  of  H  was  ten  times  larger;  this  suggested, 
rather  obviously,  that  H  was  still  the  major  element  responsible  for  the  compen¬ 
sation  of  dangling  bonds.  Second,  the  results  of  transmission  electron  microscopic 
measurements  (TEM)  on  thin  films  of  a-Si:  H,  a-Ge:  H.  a-Si,_.,Gex:  H.  a-Si:  H:  F. 
and  a-Si,_xGex:  H:  F  confirmed  earlier  indications  that  the  morphology  on  a  scale 
of  50-100  A  was  quite  likely  different.  Indeed,  the  different  structures  on  this  size 
scale  of  a-Si:  H,  a-Ge:  H  and  a-Si,_xGex:H  had  been  established  much  earlier  by  us 
[2],  and  all  subsequent  investigations  by  other  laboratories  known  to  us  have 
confirmed  our  conclusions.  The  new  element  at  this  point  is  the  finding  that  the 


alloy  a-Si,_,Ge,:  H:  F  prepared  from  fluorine-containing  gases  appears  less  defected 
than  the  corresponding  a-Si^Ge,:  H  prepared  from  the  hydride.  Our  tentative 
model,  is  that  the  presence  of  F,  or  HF,  or  other  F-containing  molecules  in  the 
preparation  plasma  can  lead  to  etching  of  the  growth  surface  of  the  film,  so  that 
weakly-coordinated  surface  entities  are  removed,  and  a  denser,  less  defected  alloy 
structure  is  produced. 

This  research  track  postulates  that  a  major  contributory  element  to  the  infe¬ 
rior  properties  of  Si^Ge,  alloys  to  those  of  Si  itself  lies  in  the  structure  on  a  100 
A  scale,  and  that  explanations  based  on  a  different  electronic  band  structure  of  a 
homogeneous  disordered  alloy  will  be  found  to  be  insufficient.  As  a  result,  we  are 
dividing  our  present  efforts  between  examinations  directed  at  investigating  the 
comparative  gap-densities-of-states  of  the  hydrogenated  and  hydro-fluorinated 
alloys  and  examinations  directed  toward  elucidation  of  the  microstructure.  In  the 
former  category,  we  include  (1)  determination  of  the  gap  density-of-states  near 
the  Fermi  level  by  the  space-charge-limited  current  method,  (2)  determination  of 
the  gap  density-of-states  by  interpretation  of  the  low  photon  energy  absorption 
spectrum  from  photoconductivity  spectra,  and  (3)  determination  by  the  time-of- 
flight  technique  of  dispersive  transport  and  deep  trapping  time.  In  the  latter 
category  we  are  investigating  (1)  the  gas  evolution  spectrum  of  a-Si,_,Ge,:  H:  F 
and  its  comparison  with  earlier  spectra  from  a-Si:  H,  a-Ge:  H  and  a-Si^Ge,:  H, 
since  different  scale  microstructure  has  usually  been  found  in  the  past  to  correlate 
with  different  photoelectronic  properties,  and  (2),  in  collaboration  with  the  group 
of  Professor  R.  Norberg  of  Washington  University  in  St.  Louis,  the  nuclear  qua¬ 
druple  resonance  of  deuterium  in  the  different  alloys  [7].  This  new  and  powerful 
technique  (for  these  alloys)  has  the  potential  for  examining  the  local  surroundings 
and  the  dynamics  of  motion  of  D  and  D2,  and  therefore  for  contributing 
significantly  to  our  knowledge  of  the  microstructure. 
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1.9  Investigation  of  Neutron  Transmutation  Doping  of  a-Si:  H.  J.H.  Bur¬ 
nett,  Z.L.  Sun,  K.D.  Mackenzie,  Y.M.  Li,  J.R.  Eggert,  and  W.  Paul,  Contract 
N00014-84-K-0465;  Research  Unit  2. 

Neutron  transmutation  doping  (NTD)  has  been  demonstrated  to  be  effective 
in  achieving  accurate,  uniform  phosphorus  doping  in  crystalline  silicon.  The  ques¬ 
tion  arises  whether  the  technique  can  be  effective  for  phosphorus-doping  hydro¬ 
genated  amorphous  silicon  and  offer  advantages  over  gas  phase  and  ion  implanta¬ 
tion  doping.  Hamanaka  et  al.  [l]  have  reported  that  NTD  is  effective  in  this  way, 
however,  since  the  results  seemed  to  us  to  be  inconclusive,  we  have  undertaken  a 
re-examination  of  this  possibility. 

Since  the  neutron  irradiation  causes  extensive  damage,  very  careful  anneal 
studies  of  the  neutron-irradiated  films  must  be  carried  out.  These  were  done  in 
parallel  with  anneal  studies  of  undoped  and  unirradiated,  and  of  light,  intermedi¬ 
ate  and  heavy  gas-phase  P  doping.  The  conductivity,  conductivity  activation 
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energy,  and  photoconductivity  of  the  different  films  where  studied  after  anneal 
times  up  to  30  hours  at  230°  C,  which  was  also  the  film  deposition  temperature. 

Based  on  a  detailed  comparison  of  the  long-time  anneal-stable  transport  and 
photo-transport,  we  have  concluded  that  the  NTD  technique  is  unsuitable  for  use 
in  a-Si:  H  technology.  At  short  times  there  do  exist  large  shifts  in  the  Fermi  level 
which  simulate  the  effect  of  considerable  P-doping,  but  these  are  related  to  defects 
and  are  not  stable.  This  study  also  showed  that  even  in  gas-phase  P-doped  a-Si:H 
there  exists  much  annealable  damage,  which  calls  into  question  the  present  view 
of  “doping  efficiency”  for  a-Si:  H.  The  conclusions  of  this  study  were  presented  at 
the  Rome  Conference  on  Amorphous  Semiconductors  [2].  There  remain  to  be 
done  other  measurements  to  consolidate  the  statistics  of  our  results  before  a  con¬ 
clusive  report  is  published. 
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Proc.  of  the  Rome  Conference  on  Amorphous  Semiconductors,  1985,  to  be 
published  in  J.  Non-Cryst.  Solids. 

1.10  Metastable  Semiconducting  Alloys  of  Ge^Sn*.  S.  Lee,  J.H.  Burnett,  and 
W.  Paul,  Contract  N00014-84-K-0465;  Research  Unit  2. 

Films  of  amorphous  Ge^Sn,  were  prepared  by  r.f.  sputtering  of  composite 
crystalline  targets  for  nine  values  of  x  between  0  and  0.4.  The  substrate  tem¬ 
perature  was  kept  low  by  circulating  water  at  room  temperature.  The  composi¬ 
tion  of  the  amorphous  films,  all  about  10  microns  thick,  was  determined  by  elec¬ 
tron  microprobe. 


At  first  the  films  were  removed  from  the  substrate  by  simply  scraping  them 
off.  Subsequent  runs  in  a  differential  scanning  calorimeter  showed  one  peak  only 
(see  below).  A  second  run  in  the  same  sample  indicated  a  feature  at  232  °  C 
clearly  identified  with  the  melting  of  £-Sn.  It  was  therefore  inferred  that  at  least 
some  of  the  material  scraped  off  the  substrate  was,  in  fact,  microcrystalline 
Ge[_jSn, .  Other  methods  of  preparing  a  free-standing  film  were  sought. 

The  technique  adopted  involved  depositing  the  Ge-Sn  alloy  onto  polyimide 
film  coating  a  microscope  glass  slide,  and  later  floating  the  film  free  by  dissolving 
the  polyimide  in  methylene  chloride.  The  largest  pieces  of  film  so  produced  were 
about  5  mms  in  diameter,  satisfactorily  large  enough  for  electrical  and  optical 
measurements. 

All  of  these  films  were  annealed  in  a  differential  scanning  calorimeter  (DSC) 
to  approximately  550  "C.  For  all  values  of  x  >  0,  two  peaks  were  observed  with 
temperature  increase:  one  inferred  to  be  caused  by  crystallization,  and  the  second 
by  phase  separation.  A  second  run  on  each  film,  increasing  the  temperature,  gave 
a  peak  at  232  0  C,  which  clearly  corresponded  to  the  melting  temperature  of  /J-Sn. 

The  temperature  of  the  first  (crystallization)  peak  varied  linearly  with  x  from 
a  value  of  500 0  C  for  x  =  0  to  an  (extrapolated)  value  of  25  0  C  for  x  =  0.5.  The 
temperature  intervals  between  the  crystallization  and  phase  separation  peaks  in 
the  DSC  did  not  depend  on  x,  but  appeared  to  vary  with  the  heating  rate  in  the 
calorimeter.  More  work  on  this  dependence  is  required  and  is  planned. 

The  actual  interval  varied  between  20  °C  and  100°  C.  From  this  result  it 
seems  likely  that  it  will  be  possible  to  anneal  the  film  to  an  intermediate  tempera¬ 
ture,  where  phase  separation  has  not  proceeded  very  far.  This  is  precisely  the 
regime  sought  by  our  experiments. 
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1.11  X-Ray  Scattering  from  Liquid  and  Solid  Surfaces.  A.  Braslau,  D. 
Schwartz,  B.  Ocko,  I.  Tidswell,  G.  Swisslow.  and  P.S.  Pershan,  Contracts 
N00014-84-K-0465,  NSF  DMR-82-12189  and  NSF  Travel  Grant  INT-83- 
11841;  Research  Unit  3. 

We  have  been  able  to  improve  the  operation  of  the  liquid  surface  spectrome¬ 
ter  that  is  installed  on  the  Harvard  MRL  Rotating- Anode  Facility  such  that  for 
some  experiments  the  signal-to-background  is  within  one  order-of-magnitude  of 
the  signal-to-background  that  has  been  obtainable  at  the  German  synchrotron 
facility  HASYLAB.  This  was  achieved  by  repositioning  certain  slits,  relative  to 
the  monochromator,  in  such  a  way  as  to  make  the  mechanical  alignment  simpler. 
This  had  the  effect  that  one  could  use  a  wider  vertical  slit  (coarser  transverse  out 
of  plane  resolution)  and  thereby  increase  the  incident  intensity.  Although  this  is  a 
considerable  improvement  for  study  of  surface  roughness  on  simple  liquids  it  is 
not  appropriate  for  study  of  the  smectic-like  order  at  liquid  crystal  surfaces  since 
the  latter  require  higher  resolution,  and  therefore  narrower  slits. 


As  a  result  we  have  been  able  to  do  much  more  in  the  way  of  sample  control 
and  development  then  we  could  when  we  were  limited  to  only  a  limited  number  of 
synchrotron  measurement  periods  in  any  given  year.  For  example,  at  the  close  of 
the  last  report  period  we  had  made  preliminary  synchrotron  measurements  on  one 
sample  of  the  AOT-oil-water  microemulsion  system  that  demonstrated  some 
structure  along  the  normal  to  the  surface  [l]. 

We  have  now  repeated  that  measurement  with  a  new  sample  on  the  Harvard 
facility  and  also  looked  at  a  second  sample  of  different  concentration.  On  the 
basis  of  both  the  repeatability  and  the  systematics  we  now  believe  that  the 
observed  structure  is  not  any  sort  of  artifact  and  in  the  next  report  period  we 
plan  to  perform  a  more  complete  set  of  measurements  and  to  also  develop 
mathematical  models  for  the  surface  properties. 

We  are  continuing  our  studies  of  the  surface  roughness  of  solid  surfaces.  We 
have  been  able  to  identify  a  number  of  empirical  phenomena  that  are  directly 
related  to  the  surface  roughness  and  we  are  in  the  process  of  trying  to  develop 
theoretical  models  to  explain  the  data.  In  collaboration  with  the  Harvard  MRL 
program  (Prof.  G.  Whitesides  of  Chemistry)  we  have  also  started  to  investigate 
the  effects  of  organic  monolayer  coatings  on  selected  solid  substrates.  These  coat¬ 
ings  have  the  potential  advantage  that  chemical  techniques  can  be  used  to 
engineer  coatings  with  diverse  properties.  People  working  on  these  types  of  coat¬ 
ings  have  been  hampered  by  the  lack  of  good  diagnostic  methods  to  evaluate  the 
chemically  produced  surfaces.  We  believe  that  x-ray  reflectivity  is  an  ideal 
method  for  this  purpose. 

Finally,  we  have  designed  a  surface  spectrometer  that  will  be  constructed  on 
beam  line  22B  on  the  x-ray  ring  at  the  National  Synchrotron  Light  Source  at 
Brookhaven  National  Laboratories.  The  Harvard  contribution  to  capital  equip¬ 
ment  costs  for  this  spectrometer  are  being  supported  bv. the  National  Science 


Foundation  through  the  Harvard  Materials  Research  Laboratory.  Travel  costs 
incurred  by  Prof.  Pershan  in  developing  the  spectrometer  are  partially  supported 
by  the  JSEP  program  and  the  NSF  under  a  separate  grant  to  Prof.  Pershan.  We 
anticipate  that  in  future  years  all  of  the  liquid  surface  synchrotron  studies  sup¬ 
ported  by  the  JSEP  grant  will  be  done  using  this  spectrometer  and  that  use  of  the 
German  facility  HASYLAB  will  be  phased  out. 

Reference: 

1.  M.  Kotlarchyk,  S-H.  Chen,  J.S.  Huang,  and  M.W.  Kim,  Phys.  Rev.  Lett.  53, 

941  (1984). 

1.12  Light  Scattering  from  Thin  Smectic  Films.  S.  Amador,  Contracts 

N00014-84-K-0465  and  NSF  DMR-82-12189;  Research  Linit  3. 

Data  collection  on  the  light  scattering  experiment  is  at  the  end  of  the  learn¬ 
ing  curve  and  yielding  results.  Problems  in  sample  purity  and  data  analysis  and 
taking  have  been  resolved,  and  data  on  DOBAMP  is  available  for  temperatures 
deep  in  the  SmC  regime.  Films  are  now  stable  into  the  transition  regions,  so  that 
critical  measurements  can  be  made  at  both  the  SmC-SmA  and  the  SmI-SmC  tran¬ 
sitions.  It  was  decided  to  do  the  SmI-SmC  transition  first  because  of  speculation 
that  surface  SmI  phase  might  form  before  the  actual  bulk  transition.  New 
interest  in  this  technique  came  about  when  synchrotron  x-ray  measurements  [l] 
revealed  the  presence  of  a  hexatic  phase  on  the  surface  of  a  free-standing  film  of 
70.7  at  temperatures  where  the  interior  has  the  bulk  Sm-C  phase.  Since  both 
these  phases  are  birefringent,  the  possibility  exists  of  measuring  pretransitional 
effects  in  this  surface  phase  using  light  scattering. 

Several  problems  were  overcome  in  the  course  of  these  measurements.  The 
laser,  a  Spectra  Physics  helium  neon  laser,  experienced  problems  due  to  aging 
vacuum  tubes  and  a  bad  plasma  tube.  These  maintenance  problems  have  been 
remedied. 
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In  the  past  year,  the  light  scattering  experiment  has  been  fully  automated  so  I 

I 

that  measurements  can  be  made  through  semi-autonomous  scans.  Data  from  the 
experiment  can  now  be  transferred  from  the  dedicated  IBM  personal  computer  to 
the  group’s  more  powerful  PDP 11/73  for  data  analysis  and  viewing.  This  necessi¬ 
tated  bringing  up  I<ERMIT,  a  communications  program,  on  both  computers,  and 
writing  macros  to  enable  use  of  standard  plotting  and  fitting  routines.  It  also 
involved  the  purchase  of  an  additional  communications  board  for  the  IBM  to 
enable  data  transfer.  In  addition,  the  minicomputer  was  a  recent  acquisition  and 
required  an  investment  of  time  in  learning  to  use  it  and  its  unfamiliar  operating 
system.  Programs  to  compensate  for  experimental  errors  are  almost  complete 
also. 

Reference: 

E.B.  Sirota,  P.S.  Pershan,  S.  Amador,  and  L.B.  Sorensen,  “X-Ray  Observa¬ 
tion  of  a  Single  Hexatic  Layer  on  the  Surface  of  a  Smectic-C  Liquid  Crystal 
Film,”  to  be  submitted  to  Phys.  Rev.  Lett.. 
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E.B.  Sirota,  P.S.  Pershan,  S.  Amador,  and  L.B.  Sorensen,  “X-Ray  Observa¬ 
tion  of  a  Single  Hexatic  Layer  on  the  Surface  of  a  Smectic-C  Liquid  Crystal 
Film,”  to  be  published. 


b.  Papers  Published  in  Refereed  Journals 

A.  Braslau,  M.  Deutsch,  P.S.  Pershan,  A.H.  Weiss,  J.  Als-Nielsen,  and  J. 
Bohr,  “Surface  Roughness  of  Water  Measured  by  X-Ray  Reflectivity,”  Phys. 
Rev  Lett  54,  114  (1985). 


1.13  Chaos  in  Driven  Josephson  Junctions.  Qing  Hu,  M.  Iansiti,  J.U.  Free, 

and  M.  Tinkham,  Contracts  N00014-83-K-0383  and  N00014-84-K-0465; 

Research  Unit  4. 

When  a  nonlinear  device  such  as  a  Josephson  junction  is  driven  periodically 
at  a  sufficient  amplitude,  the  response  may  cease  to  be  periodic,  and  acquire  a 
noisy,  aperiodic  time-dependence,  referred  to  as  “chaos.”  Since  Josephson  mixers 
must  be  driven  into  the  strongly  nonlinear  regime  to  be  efficient,  it  is  a  matter  of 
practical  concern  that  other  parameters  be  chosen  so  as  to  avoid  stumbling  into  a 
chaotic  regime,  with  its  extremely  poor  signal-to-noise  ratio.  To  be  able  to  avoid 
this  problem  in  design,  one  must  understand  the  numerous  routes  to  chaos  which 
a  given  system  may  display,  at  least  well  enough  to  avoid  them.  On  the  other 
hand,  the  study  of  these  chaotic  phenomena  is  currently  a  forefront  area  of 
research  because  of  the  insight  offered  by  model  systems  into  the  generation  of 
random  noise  from  a  noiseless  drive  in  deterministic  systems. 

During  the  current  period,  we  have  extended  [l]  the  preliminary  observations 
and  simulations  reported  earlier.  The  junctions  studied  were  Nb-aSi-Nb  SNAP 
junctions,  fabricated  by  the  Sperry  Lab.  in  Sudbury  before  it  closed,  and  supplied 
to  us.by  Dr.  L.N.  Smith.  These  junctions  are  robust  enough  to  be  studied  over  a 
long  period  of  time  without  change  of  parameter  values,  and  of  low  enough  resis¬ 
tance  relative  to  the  antenna  impedance  that  they  can  be  modeled  as  current- 
driven  at  submillimeter  wave  frequencies.  (Voltage  driven  high  R  junctions  do 
not  show  chaos,  but  they  suffer  from  noise-rounded  characteristics.)  Since  the 
plasma  frequency  of  these  junctions  is  about  400  GHz,  our  limited  FIR  (far- 
infrared)  laser  power  is  most  effective  in  driving  them  into  the  chaotic  regime 
when  a  laser  line  near  that  frequency  is  used. 

We  have  compared  our  I-V  curve  data  taken  at  a  range  of  FIR  laser  power 
levels  with  detailed  simulations,  based  on  the  RSJ  model  extended  to  take  account 
of  the  marked  change  in  quasiparticle  conductance  at  the  energy-gap  voltage.  In 


this  way,  we  have  identified  two  distinct  chaotic  regimes.  In  the  more  common  of 
the  two,  the  motion  is  characterized  by  a  “strange  attractor,”  i.e.,  even  in  the 
absence  of  any  noise  input,  the  driven  deterministic  motion  is  not  periodic,  but 
shows  a  noisy  chaotic  variation  from  cycle  to  cycle.  However,  we  have  also 
identified  chaotic  motion  in  a  parameter  range  where  the  noiseless  equations  of 
motion  lead  to  one  of  two  strictly  periodic  solutions,  corresponding  to  the  zeroth 
or  the  2/3  Josephson  step,  respectively.  The  range  of  initial  conditions  leading  to 
each  of  these  periodic  motions  is  described  as  the  “basin  of  attraction”  for  that 
motion.  In  this  case,  these  basins  are  extremely  finely  intertwined,  being 
separated  by  a  basin  boundary  of  fractal  dimension  larger  than  1  and  as  high  as 
1.85.  As  pointed  out  by  Grebogi,  McDonald,  Ott,  and  Yorke  [2],  this  fractal 
dimension  is  a  measure  of  how  sensitive  these  periodic  solutions  are  to  being  inter¬ 
rupted  by  small  amounts  of  noise  input,  since  as  the  dimension  approaches  2,  the 
boundary  fills  the  entire  phase  plane.  Our  experimental  data  confirm  this  qualita¬ 
tive  relationship,  the  dependence  on  dc  bias  current  of  the  measured  low-frequency 
noise  power  being  very  similar  to  the  dependence  on  »dc  of  the  fractal  dimension. 
To  the  best  of  our  knowledge,  this  is  the  first  direct  confirmation  of  this  plausible 
theoretical  prediction.  Moreover,  our  simulations  show  for  the  first  time  that  the 
amount  of  noise  required  to  convert  the  periodic  solution  to  a  chaotic  one,  switch¬ 
ing  intermittently  between  the  two  steps,  is  small  enough  to  be  supplied  by  the 
electronic  shot  noise  of  the  quasiparticle  current  in  the  junction,  even  if  the  device 
were  free  of  all  thermal  or  extrinsic  noise.  In  addition,  the  simulations  show  that 
the  absolute  noise  levels  resulting  from  these  two  very  different  types  of  chaos 
should  be  similar,  as  is  observed  experimentally. 

Because  the  experimental  noise  observations  are  made  at  audio  frequencies 
109  times  lower  than  the  FIR  drive  frequency  at  which  the  chaotic  motion  occurs, 
the  frequency  spectrum  of  the  noise  is  critical  to  a  quantitative  comparison  of 


simulation  and  experiment.  Our  preliminary  low-frequency  noise-spectrum  data 
show  inverse-frequency  dependences  ranging  from  f~l  to  f~*  depending  on  bias 
conditions.  A  systematic  study  of  these  dependences  is  now  underway. 

References: 

1.  M.  Iansiti,  Qing  Hu,  R.M.  Westervelt,  and  M.  Tinkham,  Phys.  Rev.  Lett.  55, 
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M.  Iansiti,  Qing  Hu,  R.M.  Westervelt,  and  M.  Tinkham,  “Noise  and  Chaos  in 
a  Fractal  Basic  Boundary  Regime  of  a  Josephson  Junction,”  Phys.  Rev.  Lett. 
55,  746-749  (1985). 
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Conferences 

M.  Tinkham,  “Small-Area  Josephson  Devices  for  Submillimeter 
Wavelengths,”  Invited  Paper  at  US-Japan  Workshop  on  Josephson  Electron¬ 
ics,  Kauai,  Hawaii,  June,  1985. 


1.14  Nonlinear  Dynamics  of  Small  Neural  Networks  with  Finite  Fre¬ 
quency  Response.  K.  Babcock  and  R.M.  Westervelt,  Contract  N00014-84- 
K-0465;  Research  Unit  5. 

We  are  investigating  the  nonlinear  dynamics  of  networks  of  electrical  thres¬ 
hold  switching  elements  which  mimic  biological  neural  networks.  These  show  con- 
siderable-promise  for  rapid  pattern  recognition  and  optimization  problems;  how¬ 
ever,  their  dynamics  resemble  those  of  a  spin  glass  and  are  only  poorly  under¬ 
stood.  Most  work  to  date  has  focussed  on  algorithms  to  produce  multiple  steady 
states  in  simulations  of  large  networks.  Serious  performance  problems  may  occur 
in  real  networks  where  the  switching  elements  have  finite  frequency  response, 
including  the  appearance  of  unwanted  collective  modes  produced  by  phase  shifts 
in  real  elements,  and  exceptional  sensitivity  to  noise  produced  by  complex  basin 
boundaries  [l] .  The  significance  of  these  problems  can  be  appreciated  by  the 
difficulty  of  operating  a  conventional  digital  computer  without  a  synchronizing 
clock. 

We  have  investigated  the  dynamics  of  relative  simple  networks,  using 
mathematical  analysis,  computer  simulations,  and  experiments  on  circuits  com¬ 
posed  of  integrated  circuit  amplifiers  [2].  Our  goal  has  been  to  understand  the 
nonlinear  dynamics  of  small  systems  with  2  to  4  nodes  (and  4  to  16  degrees  of 
freedom)  in  some  detail,  in  order  to  identify  simple  characteristic  phenomena 


which  can  be  used  to  analyze  the  behavior  of  larger  networks  with  greater  poten¬ 
tial  computational  power,  but  hopelessly  complex  dynamics. 

In  a  network  with  two  switching  elements  and  two  steady  states,  we  find 
rapid  convergence  to  the  steady  states,  and  a  smooth  boundary  separating  their 
basins  of  attraction  only  if  the  frequency  response  of  the  elements  is  assumed  to 
be  infinitely  fast.  When  a  more  realistic  model  including  phase  shifts  in  each  ele¬ 
ment  is  used,  the  response  of  the  network  is  much  more  complex.  The  transient 
response  displays  damped  ringing  even  when  the  interconnection  matrix  T,-,  is 
symmetric,  and  the  response  to  periodic  excitation  becomes  irregular  and  chaotic 
in  certain  cases  when  the  drive  frequency  is  near  the  high  frequency  cut  off  of  the 
elements.  Under  these  circumstances  one  also  expects  that  the  boundaries  of  the 
basins  of  attraction  for  the  two  steady  states  are  highly  folded  fractal  sets  which 
cause  extraordinary  sensitivity  to  the  choice  of  initial  conditions;  this  issue  is 
currently  under  investigation.  A  small-signal  analysis  defines  potentially  trouble¬ 
some  combinations  of  operating  parameters  and  provides  a  means  for  avoiding 
irregular  response. 
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S.W.  Teitsworth  and  R.M.  Westervelt,  “Subharmonic  and  Chaotic  Response 
of  Periodically-Driven  Extrinsic  Ge  Photoconductors,  submitted  to  Phys 
Rev.  Lett. 
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in  the  Damped  Driven  Pendulum,”  in  Proc.  17th  Int.  Conf.  on  Low  Tempera¬ 
ture  Physics,  Karlsruhe,  W.  Germany,  U.  Eckern,  A.  Schmid,  W.  Weber  and 
H.  Wuhl,  eds.,  (North-Holland,  Amsterdam,  1984),  p.  1139. 
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E.G.  Gwinn  and  R.M.  Westervelt,  “Intermittent  Chaos  and  Low  Frequency 
Noise  in  the  Driven  Damped  Pendulum,”  Phys.  Rev.  Lett.  54,  1613  (1985). 


Books  (and  sections  thereof)  Submitted  for  Publication 
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linear  Dynamics  (World  Scientific  Publishers,  Singapore,  1986),  in  press. 


b  y 


-  24  - 


g.  Invited  Presentations  at  Topical  or  Scientific/Technical  Society 
Conferences 

R.M.  Westervelt,  “Chaos  in  Semiconductors  and  in  Josephson  Junctions,” 
Universite'  de  Sherbrooke,  Regional  School  of  Physics,  29-31  October,  1984. 
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n.l  Experimental  Studies  of  Photon  Echo  Pulse  Compression.  Y.S.  Bai 
and  T.W.  Mossberg,  Contracts  N00014-84-K-0465  and  ARO  DAAG29-83-K- 
0040;  Research  Unit  6. 

Following  up  on  our  previous  theoretical  work  [l],  we  have  performed  an 
experimental  study  of  the  compression  of  frequency-chirped  optical  pulses  using 
photon-echo-type  effects  [2].  We  modified  our  cw  ring  dye  laser  so  that  we  could 
generate  long  (hundreds  of  nanoseconds)  pulses  possessing  a  nearly  linear  fre¬ 
quency  chirp.  Photon  echoes  were  generated  by  passing  two  frequency-chirped 
pulses  through  an  oven  containing  vapor  phase  atomic  Yb.  When  the  chirp  rate 
of  the  second  laser  pulse  was  adjusted  to  be  twice  that  of  the  first,  the  duration  of 
the  echo  signal  was  found  to  be  given  by  the  inverse  of  the  total  chirp  bandwidth. 
Echoes  having  durations  up  to  30  times  shorter  than  the  first  excitation  pulse 
were  observed.  Moreover,  the  echo  signal  was  found  to  contain  up  to  25^0  of  the 
energy  contained  in  the  first  laser  excitation  pulse.  The  present  means  of 
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compressing  optical  pulses  is  unique  in  that  long  input  pulses  can  be  compressed 
as  well  as  short  input  pulses.  Other  techniques  of  optical  pulse  compression 
become  difficult  to  employ  when  the  input  pulses  are  long. 
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II. 2  Broad-Bandwidth  Transient  Four-Wave  Mixing.  J.E.  Golub  and 
T.W.  Mossberg,  Contracts  N00CH.4-84-K-0465  and  NSF  PHY-82-07080; 
Research  Unit  6. 

It  was  recently  suggested  [l]  that  transient  nonlinear  optical  experiments 
employing  broad-bandwidth  laser  pulses  may  be  used  to  provide  ultra-high  tem¬ 
poral  resolution  in  optical  studies  of  material  systems.  In  these  experiments,  •• 
mixing  signal  is  measured  as  a  function  of  an  optical  delay  introduced  into  the 
path  of  one  of  the  input  beams  (at  least  two  of  which  must  be  optically  divided 
portions  of  the  same  laser  output  pulse).  The  mixing  signal  has  been  calculated  to 
decay  on  the  time  scale  of  the  longer  of  the  laser  coherence  time  or  the  transverse 
relaxation  time  of  the  material  employed.  The  overall  laser  pulse  duration  is  not 
important  unless  the  optical  delay  becomes  comparable  to  it.  As  a  result,  meas¬ 
urements  of  material  relaxation  rates  should  be  limited  only  by  the  coherence  time 
of  the  laser  pulse  employed. 

We  have  conducted  an  experiment  [2]  in  a  simple  collisionally  perturbed  gas- 
phase  system  designed  to  test  the  ideas  described  above.  A  10-nanosecond  dye 
laser  pulse  with  a  1-10  nm  bandwidth  was  optically  split  into  three  parts  and  used 
to  generate  a  four-wave  mixing  signal  in  atomic  sodium  vapor.  By  using  high 
argon  perturber  pressures  (1-4  atmospheres)  we  were  able  to  introduce  collisional 


decay  times  of  10-30  picoseconds.  We  found  that  the  four-wave  mixing  signal 
decayed  exponentially  as  a  function  of  optical  decay,  and  that  the  picosecond 
decay  constant  agreed  beautifully  with  the  collisional  decay  rate  expected  on  the 
basis  of  traditional  collision-broadening- type  experiments.  Our  experiment,  per¬ 
formed  in  a  well  understood  material  system,  is  essentially  the  first  to  show  that 
broad-bandwidth  transient  mixing  techniques  can  indeed  provide  reliable  measure 
ments  of  ultra-high  speed  relaxation  processes. 
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ing,”  submitted  to  the  J  Opt.  Soe  B. 

II.3  Interaction  of  Atoms  with  Amplitude-  and  Phase-Controlled  Opti¬ 
cal  Radiation.  Y.S.  Bai,  A.G.  Yodh.  and  T.W.  \.  >s>herg.  Contracts 
N000 1 4-84-K-04 65  and  ONR  N00014-85-K-O72 1:  Research  Unit  6. 

Although  the  behavior  of  two-level  atoms  exposed  to  nearly  resonant  laser 
light  can  be  simply  calculated  in  the  semiclassical  approximation  and  is  generally 
assumed  to  be  well  understood,  no  experiments  have  been  performed  to  date  to 
actually  observe  the  response  of  such  atoms  under  conditions  in  which  the  atom- 
laser  relative  phase  is  important.  In  typical  steady-state  experiments,  of  which 
there  are  many,  the  observed  signal  depends  on  the  ensemble  averaged  behavior  o 
atoms  with  a  wide  range  of  phases  relative  to  the  applied  field.  In  transient 
experiments,  the  laser’s  phase  is  generally  held  constant. 

We  have  begun  a  series  of  experiments  designed  to  probe  the  behavior  of 
atoms  responding  to  laser  fields  of  well-specified  phase  [  1  ] .  These  experiments  are 
carried  out  in  dilute  atomic  beams  of  atomic  Yb.  When  excited  by  circularly 


polarized  555.6  nm  light,  this  atom  constitutes  an  ideal  two-level  system.  In  our 
experiments,  atoms  initially  in  the  ground  state  are  excited  by  a  gated  cw  laser 
field  having  a  controllable  phase.  We  monitor  the  atoms  behavior  by  observing 
fluorescence  emitted  at  right  angles  to  both  the  atomic  and  laser  beams. 

In  initial  experiments  [2],  we  have  studied  the  behavior  of  atoms  that  are 
first  exposed  to  a  pulse  of  area  ir/2  (i.e.,  they  are  placed  in  an  equal  admixture 
superposition  of  ground  and  excited  states)  and  then  exposed  to  a  laser  field 
shifted  in  phase  by  90  degrees  with  respect  to  the  initial  field.  We  have  pointed 
out  that  this  excitation  process  places  the  combined  atom-laser  field  system  into  a 
stationary  state.  In  resonance  fluorescence  theory,  these  states  are  called  “dressed 
states.”  We  have  also  pointed  out  that  the  resonance  fluorescence  spectrum 
characteristic  of  atoms  excited  to  particular  dressed  states  will  be  qualitatively 
different  from  that  normally  observed  in  steady-state  experiments.  We  have 
demonstrated  that  atoms  in  dressed  states  do  not  exchange  energy  with  the  exci¬ 
tation  field. 

These  experiments  probe  basic  aspects  of  the  interaction  of  simple  atomic 
systems  with  laser  radiation. 
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H. 4  Collisional  Relaxation  of  the  3PX  State  of  Atomic  Yb.  A.G.  Yodh  and 

T.W.  Mossberg,  Contracts  N00014-84-K-0465  and  NSF  PHY-82-07080; 
Research  Unit  6. 

We  have  employed  [l]  coherent  transient  techniques  to  study  the  relaxation 
of  the  3Pj  state  of  atomic  Yb  collisionally  perturbed  by  helium  and  argon. 
Specifically,  our  measurements  are  sensitive  to  velocity  and  phase  thermalization 
of  the  rank  one  (orientation)  and  rank  two  (alignment)  multipole  moments  of  the 
3P,  Zeeman  manifold.  We  provided  the  first  measurement  of  the  total  quantum 
mechanical  cross-section  for  destruction  of  excited-state  Zeeman  coherences  and 
the  first  complete  mapping  of  the  collisional  velocity-change  spectrum  characteris¬ 
tic  of  the  collisions. 
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A32,  844  (1985). 

II. 5  Mixed  Binary  Multiplication  Using  Real-Time  Optical  Convolution. 
W.R.  Babbit  and  T.W.  Mossberg,  Contracts  N00014-84-K-0465  and  ARO 
DAAG29-83-K-0040;  Research  Unit  6. 

Using  the  method  of  real-time  optical  convolution  developed  here  previously 
[l],  we  have  shown  [2]  that  it  is  possible  to  perform  nanosecond-time-scale  mixed 
binary  multiplication  of  serially  encoded  optical  data.  In  this  process,  t  wo  or 
three  numbers  are  digitally  encoded  as  an  optical  pulse  train  and  made  to  interaci 
with  a  resonant  inhomogeneously  broadened  optical  absorber.  The  input  pulses 
modify  the  absorption  profile  of  the  material  to  represent  the  product  of  their 
Fourier  transforms.  Subsequent  coherent  free  decay  signals  emitted  by  the  sample 
effectively  reverse  transform  the  stored  information  and  hence  evolve  in  time  as 
the  convolution  or  cross-correlation  of  the  input  pulses.  Mixed  binary  multiplica¬ 
tion  occurs  when  two  or  more  binary  encoded  signals  are  convolved.  Our 


experiments  indicate  that,  in  weakly  relaxing  materials,  the  mixed  binary  output 
signal  obtained  with  the  present  technique  should  be  about  0.5%  as  intense  as  the 
input  signals.  Inaccuracies  in  the  output  signals  were  investigated  and  ascribed  to 
correctable  imperfections  in  our  experimental  system.  Harvard  has  applied  for  a 
patent  on  this  process. 
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J.E.  Golub  and  T.W.  Mossberg,  “Terahertz  Quantum  Beat  Spectroscopy,” 
presented  at  the  1985  Gordon  Conference  on  Lasers  and  Nonlinear  Optics, 
August  1,  1985). 

A.G.  Yodh,  Y.S.  Bai,  and  T.W.  Mossberg,  “Studies  of  Two-Level  Atoms 
Identically  Prepared  by  Phase  and  Amplitude  Controlled  Excitation,” 
presented  at  the  1985  Gordon  Conference  on  Atomic  Physics,  July  1-5,  1985. 

A.G.  Yodh,  J.E.  Golub,  and  T.W.  Mossberg,  “The  Collisional  Relaxation  of 
Excited-State  Zeeman  Coherences  in  Atomic  Ytterbium  Vapor,”  presented  at 
the  XTV  International  Conference  on  the  Physics  of  Electronic  and  Atomic 
Collisions,  Palo  Alto,  CA,  July  23  -  August  1,  1985. 
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Y.S.  Bai  and  T.W.  Mossberg,  “Photon  Echo  Pulse  Compression,”  presented 
at  the  1985  Conference  on  Lasers  and  Electro-Optics  (CLEO),  May  22,  1985. 

W.R.  Babbitt,  Y.S.  Bai,  N.W.  Carlson,  and  T.W.  Mossberg,  “Spectral  Holog¬ 
raphy  and  its  Applications,”  presented  at  the  1985  Conference  on  Lasers  and 
Electro-Optics,  May  22,  1985. 

Y.S.  Bai,  A.G.  Yodh,  T.W.  Mossberg,  N.  Lu,  and  P.R.  Berman,  “Transient 
Excitation  of  Atoms  Using  Strong  Phase-Controlled  Laser  Fields,”  presented 
at  the  1985  Gordon  Conference  on  Lasers  and  Electro-Optics,  August  1, 

1985. 

Y.S.  Bai,  A.G.  Yodh,  and  T.W.  Mossberg,  “Studies  of  Two-Level  Atoms 
Identically  Prepared  by  a  Phase-  and  Amplitude-Controlled  Excitation 
Field,”  presented  at  the  Seventh  International  Conference  on  Laser  Spectros¬ 
copy,  Hawaii,  June  24-28,  1985. 

A.G.  Yodh,  J.  Golub,  and  T.W.  Mossberg,  “Colliding  Without  Relaxing:  The 
Suppression  of  Collisional  Dephasing  with  Strong  Optical  Fields,”  presented 
at  the  Seventh  International  Conference  on  Laser  Spectroscopy,  Hawaii,  June 
24-28,  1985. 


H.6  Nonlinear  Four-Wave  Mixing  in  Vapors.  N.  Bloembergen  and  Y.H. 

Zou,  Contract  N00014-84-K-0465;  Research  Unit  7. 

During  the  past  year  a  comprehensive  experimental  study  of  the  collision- 
enhanced  Hanle-type  resonances  in  Na  vapor  with  various  buffer  gases  has  been 
completed.  After  the  original  communication  [1],  referred  to  in  the  preceding 
annual  report,  the  intensity  and  width  of  the  zero  magnetic  field  resonances  in 
phase-conjugate  degenerate-frequency  four-wave  light  mixing  have  been  systemati¬ 
cally  investigated  as  a  function  of  several  parameters.  These  include  the  pressure 
of  the  Ar  and  Xe  buffer  gas,  the  partial  pressure  of  Na  vapor,  the  detuning  of  the 
light  frequency  on  both  sides  of  the  Z?t  and  D2  sodium  resonances  and  the  inten¬ 
sity  of  the  incident  light  beams.  Power  broadening  and  sat  uration  effects  have 
been  observed  and  interpreted.  It  has  been  demonstrated  theoretically  that  the 
collision-enhanced  Hanle  resonances  are  equivalent  to  the  phenomenon  of 
collision-induced  transverse  optical  pumping.  The  relationship  between  the  two 
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theoretical  descriptions  simply  corresponds  to  a  90  degree  rotation  of  the  axis  of 
spatial  quantization.  The  detuning  from  resonance  has  been  extended  well  over 
100  GHz  to  about  600  GHz.  In  this  regime  the  impact  approximation  is  not  valid 
and  the  line  shape  becomes  asymmetrical.  This  is  the  regime  of  the  quasi-static 
collisional  deformation  of  the  wave  functions.  All  these  effects  are  described  in 
detail  in  a  comprehensive  paper,  which  has  been  submitted  for  publication  [2]. 

Some  preliminary  experiments  have  been  carried  out  on  the  collision- 
enhanced  Hanle  resonances  in  Samarium  vapor.  A  new  oven  has  been  constructed 
which  permits  heating  at  800  °  C  to  obtain  sufficient  vapor  pressure  of  Sm  atoms. 
This  atom  is  of  interest  because  it  avoids  the  complications  of  the  presence  of 
hyperfine  interaction  present  in  Na.  Furthermore,  the  near-resonant  excited  state 
is  a  single  ( J  =  o)  state,  while  the  ground  state  is  a  triply  degenerate  [J  =  l) 
state  in  zero  magnetic  field.  The  systematic  variation  of  the  states  of  polarization 
in  the  simple  energy  level  scheme  is  of  special  interest.  The  Hanle  resonance  of 
the  Sm  atom  with  Ar  and  Xe  buffer  gases  has  been  observed.  A  systematic  inves¬ 
tigation  of  the  parametric  dependences  of  the  intensity  and  width  of  this  Hanle- 
type  signal  is  underway. 
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H.7  Picosecond  Laser  Interactions  with  Semiconductors.  A.M.  Malvezzi, 

C.Y.  Huang,  and  N.  Bloembergen,  Contracts  N00014-84-K-0465  and 

N00014-83-K-0030;  Research  Unit  7. 

During  this  reporting  period  the  experiments  on  the  emission  of  neutral  and 
charged  particles  from  Si  and  GaAs  were  continued.  For  heating  above  the  nor¬ 
mal  boiling  point  (or  vaporization  temperature),  a  regime  of  plasma  formation  is 
rapidly  reached.  When  positive  ions  are  extracted  by  applying  a  negative  bias  to 
the  collecting  electrode,  a  strong  negative  space-charge  region  is  rapidly 
developed,  which  leads  to  further  ionization  in  the  vapor  phase.  Further  theoreti¬ 
cal  studies  of  the  emmission  of  positive  ions  are  in  progress. 

Our  technique  of  a  picosecond  pump  heating  pulse  with  subsequent  probing 
of  the  changes  in  the  optical  complex  index  of  refraction  by  picosecond  probing 
pulses  at  several  wavelengths  with  various  time  delays  has  been  extended  to  gra¬ 
phite.  A  phase  transition  has  been  observed  at  a  critical  pump  fluence  of  about 
160  mJ/cm2  at  n  =  0.53  nm.  An  abrupt  change  in  the  reflectivity  has  been  observed 
at  several  wavelengths,  including  i  =  0.53  fim,  1.06  jim  and  1.9  /zm.  The  transition 
is  tentatively  identified  with  molten  carbon.  Since  the  reflectivity  decreases,  our 
results  would  indicate  that  the  electronic  structure  of  liquid  carbon  is  not  metal¬ 
lic,  but  rather  molecular  polymeric  in  nature.  Further  experimental  and  theoreti¬ 
cal  work  on  this  problem  is  in  progress. 

Our  pump  and  probe  technique  will  also  be  used  to  investigate  other  high 
melting  point  refractory  compounds,  including  silicon  and  tungsten  carbide. 

Several  papers  on  these  subjects  [l-4]  have  been  accepted  for  presentation  at 
the  Materials  Research  Society  Symposium  on  Beam-Solid  Interactions  and  Phase 
Transformations,  held  jn  Boston  in  December  1985.  They  will  also  be  submitted 
for  publication  in  the  Proceedings  of  this  symposium.  A  paper  describing  a  simple 
model  for  the  variation  of  the  effective  mass  wih  the  electron-hole  plasma  density 
and  temperature  in  silicon  has  been  accepted  for  publication  [5]. 
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Produced  High-Density  Plasma  in  Silicon,”  accepted  for  publication  in  IEEE 
Journal  of  Quantum  Electronics,  January  1986. 

b.  Papers  Published  in  Refereed  Journals 

N.  Bloembergen,  Y.H.  Zou,  and  L.J.  Rothberg,  “Collision-Induced  Hanle 
Resonances  of  Kilohertz  Width  in  Phase-Conjugate  Four-Wave  Light  Mix¬ 
ing,”  Phys.  Rev.  Lett.  54,  186-188  (1985). 

A.M.  Malvezzi,  H.  Kurz,  and  N.  Bloembergen,  “Nonlinear  Photoemission 
from  Picosecond  Irradiated  Silicon,”  Appl.  Phys.  A36,  143-146  (1985). 

H.  Kurz  and  N.  Bloembergen,  “Picosecond  Photon-Solid  Interaction,”  Mat. 
Res.  Soc.  Symp.  Proc.  Vol.  35,  Elsevier  Science  Publishing  Companv,  1985, 
pp.  1-13. 


J.M.  Liu,  A.M.  Malvezzi,  and  N.  Bloembergen,  “Ultrafast  Phase  Transition 
of  GaAs  Surfaces  Irradiated  by  Picosecond  Laser  Pulses,”  Mat.  Res.  Soc. 


Symp.  Proc.  Vol.  35,  Elsevier  Science  Publishing  Company,  1985,  pp.  137- 
142. 

A.M.  Malvezzi,  H.  Kurz,  and  N.  Bloembergen,  “Picosecond  Photoemission 
Studies  of  the  Laser-Induced  Phase  Transition  in  Silicon,”  Mat.  Res.  Soc. 
Symp.  Proc.  Vol.  35,  Elsevier  Science  Publishing  Company,  1985,  pp.  75-80. 

N.  Bloembergen  and  Y.H.  Zou,  “The  Relationship  Between  Collision-Assisted 
Zeeman  and  Hanle  Resonances  and  Transverse  Optical  Pumping,”  to  be  pub¬ 
lished  in  Laser  Spectroscopy  V ,  T.  Hansch,  Jr.  and  Y.R.  Shen,  eds.  (Springer, 
Heidelberg,  1985),  pp.  186-191. 


g.  Invited  Presentations  at  Topical  or  Scientific/Technical  Society 
Conferences 

N.  Bloembergen,  “Collision-Induced  Hanle  Resonances  in  Four-Wave  Light 
Mixing.”  Invited  Lecture  at  the  International  Symposium  in  honor  of  the 
memory  of  A.  Kastler,  Paris,  France,  January  1985. 

N.  Bloembergen,  “The  Prehistory  and  Future  of  Nonlinear  Optics.”  Invited 
Lecture  at  the  CLEO  meeting,  Baltimore,  Maryland,  May  1985,  Symposuim 
for  the  25th  Anniversary  of  the  Laser. 

N.  Bloembergen,  “The  Relationship  Between  Collision-Assisted  Zeeman  and 
Hanle  Resonances  and  Transverse  Optical  Pumping.”  Invited  Lecture  at  the 
seventh  International  Conference  on  Laser  Spectroscopy,  Maui,  Hawaii,  1985. 

N.  Bloembergen,  “Collision-Induced  Coherence.”  Invited  Led  ure  at  the  Sum- 
posium  on  Coherence,  Cooperation  and  Fluctuations,  in  honor  of  R.J. 
Glauber,  Cambridge,  Mass.,  October  1985. 


n.9  Multiphoton  Vibrational  Excitation  of  Molecules.  E.  Mazur,  J.  Wang, 
and  K.H.  Chen,  contract  N00014-84-K-0465;  Research  Unit  8. 

Time-resolved  Raman  spectroscopy  of  gases  at  very  low  densities  is  used  to 
probe  the  intramolecular  distribution  of  energy  in  multiphoton  excited  molecules. 
This  distribution  is  studied  as  a  function  of  time,  infrared  laser  fluence  and  molec¬ 
ular  size.  The  total  overhaul  of  the  experimental  set-up,  which  was  started  in  the 
previous  reporting  period,  has  been  completed  and  the  first  runs  with  the  new 
set-up  have  been  performed. 

The  optical  free  induction  decay  cell,  which  is  used  to  generate  short 
picosecond  infrared  pulses  of  controllable  duration,  was  replaced  and  the  plasma 
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shutter  improved.  A  mode-beat  detection  and  correction  system  will  soon  be 
added  to  the  oscillator.  In  order  to  substantially  reduce  the  consumption  of  gas,  a 
gas-recirculation  system  was  added  to  the  high  pressure  C02-laser  amplifier.  The 
unit  can  now  operate  without  any  flow  of  gas  for  a  period  up  to  a  day  without 
substantial  loss  of  gain.  This  will  also  allow  us  to  operate  with  isotopic  gas  mix¬ 
tures  if  needed,  e.g.  to  obtain  a  higher  gain  at  the  weak  lines  of  ordinary  C02  gas 
mixtures.  A  totally  new  computer-controlled  data-acquisition  system,  consisting 
of  a  microcomputer  and  a  CAMAC  (computer  automated  measurement  and  con¬ 
trol)  crate,  has  been  interfaced  to  the  experiment.  The  software  has  been  com¬ 
pletely  rewritten. 

After  the  overhaul  was  completed  around  October,  we  have  first  tested  the 
reproducibility  of  the  apparatus  by  repeating  previously  reported  and  published 
measurements  on  SF„  [lj.  The  measurements  reproduced  the  results  obtained  pre¬ 
viously,  demonstrating  the  reliability  of  the  apparatus.  As  planned  we  proceeded 
then  by  extending  our  study  to  other  molecular  systems.  For  CF2HC£  preliminary 
tests  in  the  old  apparatus  had  shown  that  five  Raman  active  modes  of  widely 
different  energy  can  be  observed  at  densities  low  enough  to  ensure  a  collisionless 
excitation.  This  molecule,  however,  absorbs  on  the  edge  of  the  C02-laser  range  at 
the  9.4  /im  R(32)  line.  At  the  maximum  fluence  we  could  achieve,  photoacoustic 
measurements  show  that  the  molecules  absorb  less  than  one  photon  per  molecule. 
At  those  fluences  and  excitation  levels,  no  measurable  change  in  the  Raman  spec¬ 
trum  could  *  e  observed. 

VVe  are  presently  carrying  out  measurements  on  C2H4F2.  This  molecule  has 
four  strong  Raman  active  modes,  ranging  in  energies  from  400  to  3000  cm-1,  which 
can  be  observed  in  our  apparatus.  In  contrast  to  CF2HC£,  however,  it  absorbs 
exactly  at  the  peak  of  the  C02-laser  range,  and  we  are  able  to  cover  the  entire 
excitation  range  from  ground  state  up  to  the  dissociation  limit.  When  the 


molecule  dissociates  a  complication  arises  because  the  dissociation  fragments 
absorb  the  ultraviolet  from  the  probing  pulse  and  subsequently  emit  a  broadband 
radiation,  which  extends  into  the  anti-Stokes  range  (this  implies  that  the  frag¬ 
ments  contain  a  considerable  amount  (>  3000  cm-1)  of  internal  energy  after  dissoci¬ 
ation).  By  monitoring  the  signal  in  between  two  Raman  lines  it  is  possible  to 
determine  exactly  at  what  infrared  laser  fiuence  dissociation  starts  to  play  a  role. 
These  measurements  are  being  carried  out  presently.  Subsequently  we  will  limit 
our  study  of  the  Raman  signals  for  this  molecule  to  the  range  of  fluences  below 
the  dissociation  threshold.  In  addition,  the  collisional  decay  of  the  excitation  of 
the  dissociation  fragments— which  was  observed  in  our  first  runs  and  which  was 
seen  to  vary  over  different  parts  of  the  spectrum-will  be  studied  in  more  detail. 
Reference: 

1.  E.  Mazur,  I.  Burak,  and  N.  Bloembergen,  “Collisionless  -Vibrational  Energy 
Redistribution  Between  Infrared  and  Raman  Active  Modes  in  SF„,”  Chem. 
Phys.  Lett.  105,  258  (1984). 

ANNUAL  REPORT  OF 

PUBLICATIONS/PATENTS/PRESENTATIONS /HONORS 

a.  Papers  Submitted  to  Refereed  Journals  (and  not  yet  published) 

D.  Korff  and  E.  Mazur,  “Beam  Combining  Considerations  in  Stimulated 
Raman,”  J.  Opt.  Soc.  America. 

b.  Papers  Published  in  Refereed  Journals 

H.  van  Houten,  E.  Mazur,  and  J.J.M.  Beenaker,  “The  Temperature  Depen¬ 
dence  of  Flow  Birefringence  in  Gases  and  the  Scalar  Factor  of  Angular 
Momentum  Polarization  in  Viscous  Flow,”  Chem.  Phys.  Lett.  113,  135 
(1985). 

E.  Mazur,  H.  van  Houten,  and  J.J.M.  Beenaker,  “A  Comparison  of  Data  on 
the  Viscomagnetic  Effect,  Flow  Birefringence  and  Depolarized  Rayleigh  Line 
Broadening,”  Physica  130A,  505  (1985). 
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E.A.  Mason  and  E.  Mazur,  “Field  Effects  on  Transport  Phenomena  in  the 
Transition  Regime,”  Physica  120A,  437  (1985). 


Honors  /Awards /P  r  izes 

E.  Mazur,  Inclusion  in  Who’s  Who  in  Frontiers  of  Science  and  Technology, 
1985. 
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in.  INFORMATION  ELECTRONICS  CONTROL 
AND  OPTIMIZATION 
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Prof.  R.W.  Brockett 

Prof.  Y.C.  Ho 
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Asst.  Prof.  R.  Suri 

Dr.  X.  Cao 
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Dr.  J.  Loncaric 

Dr.  A.  Mclvor 

Dr.  I.  Nikolic 

Dr.  R.  Rubenstein 

Dr.  K.  Wohn 


Mr.  W.  Cebuhar 
Mr.  D.  Cochran 
Mr.  J.  Dille 
Mr.  V.  Eng 
Mr.  P.  Glassman 
Mr.  W.  Gong 
Mr.  M.  Lee 
Mr.  D.  Montana 
Mr.  Y.  Wang 
Mr.  J.  Wu 
Mr.  M.  Zananis 


in.  1  Conditional  Density  Approach  to  Tracking.  R.W.  Brockett,  A. 

Mclvor,  K.  Wohn,  and  J.  Wu,  Contracts  N00014-84-K-0465,  N00014-84-K- 

0504,  DAAC29-83-K-0027  and  AFOSR-81-0054;  Research  Unit  9. 

The  problem  of  determining  shape  and  motion  parameters  from  a  sequence  of 
images  is  basic.  We  have  shown  how  to  formulate  problems  of  this  type  in  terms 
of  the  conditional  probability  density  with  the  conditioning  being  with  respect  to 
the  observations.  An  unexpected  outcome  of  this  line  inquiry  is  a  new.  and  rather 
compelling,  argument  for  the  use  of  edge  information,  both  in  conditional  density 
propagation  and  in  simplifying  the  algorithms.  We  show  that  in  a  certain  sense, 
the  complexity  of  tracking  an  object  as  well  as  possible  grows  as  the  length  of  its 
perimeter  and  not  as  its  area.  A  second  aspect  of  this  work  involves  the  develop¬ 
ment  of  statistical  models  for  the  random  aspects  of  these  problems.  Here  we 


have  concentrated  on  questions  involving  the  spatial  correlation,  or  power  spec¬ 
trum,  of  the  various  signals  involved,  including  a  stochastic  description  of 
reflectance.  This  has  led  to  an  interpretation  of  some  well  known  edge  detection 
filters  as  “prewhitening  filters”  in  the  sense  of  Bode  and  Shannon. 

Reference: 

1.  Alan  Mclvor,  Stochastically  Based  Vision  Algorithms,  Ph.D.  Thesis,  Harvard 
University,  June,  1985. 

III.2  Analysis  of  Moving  Images.  R.  Brockett,  Contracts  N00014-84-K-0504 
and  N00014-84-K-0465;  Research  Unit  9. 

The  first  issue  is  that  of  obtaining  the  optical  flow  (the  projection  of  3-D 
velocity  vector)  from  the  sequence  of  time-varying  images.  So  far  several 
researchers  have  proposed  flow  computation  methods  which  involve  tracking  a 
finite  number  of  feature  points  Which  (hopefully)  correspond  to  the  structural 
extrema  of  objects.  Since  our  motion  analysis  scheme  requires  not  only  the  flow 
itself  but  also  its  spatial  derivatives,  such  methods  are  not  adequate  for  obtaining 
accurate  flow  measurements  and  we  are  currently  investigating  two  different 
methods  as  follows: 

Significant  intensity  changes  usually  correspond  to  the  texture  boundaries 
and  physical  edges.  By  observing  the  deforming  contours  over  time,  one  can  com¬ 
pute  the  flow.  We  have  shown  that  the  flow  and  its  derivatives  up  to  second 
order  can  be  recovered  quite  reliably,  although  the  flow  up  to  the  first  order  alone 
will  suffice.  We  have  also  found  that  the  flow  may  be  much  more  accurately 
recovered  by  utilizing  the  notion  of  “temporal  smoothness”  of  flow.  According  to 
this  model  the  optical  flow  is  approximated  locally  as  the  first-order  polynomial  in 
spatio-temporal  coordinates. 

There  are  many  situations  in  which  contours  cannot  be  extracted  or  are  very 
sensitive  to  noise  in  image.  We  are  studying  a  complementary  approach  which 


relies  on  the  gray  level  changes.  Under  the  assumption  of  convected  invariance  of 
image  intensity,  temporal  change  of  intensity  can  be  measured  by  the  two 
mechanisms:  intensity  correlation  and  spatio-temporal  intensity  gradient-based. 
We  try  to  obtain  the  affine  transformation  of  local  image  deformation,  directly 
from  the  intensity  change.  Such  affine  transformation  will  determine  the  local 
surface  of  objects  and  the  global  motion  parameters  as  explained  in  the  following 
section. 

The  above  two  different  approaches  may  be  combined  in  order  to  estimate  a 
more  reliable  flow  field.  We  propose  a  multi-resolution  approach  which  (i)  con¬ 
structs  a  pyramid  of  individual  raw  image  frames,  (ii)  applies  the  individual 
methods  of  flow  recovery,  (iii)  determines  which  method  provides  the  most  accu¬ 
rate  estimates,  and  (iv)  produces  the  reliable  flow  field  by  communication  through 
the  various  levels  of  pyramid. 

Our  approach  is  based  on  the  optical  flow  as  is  viewed  as  Eulerian  description 
(opposed  to  Lagrangian).  It  has  been  shown  that  the  3-D  motion  and  object  struc¬ 
ture  (in  terms  of  surface  metric  and  curvature  tensors)  may  be  recovered  from  the 
optical  flow  field  and  its  first  and  second  derivatives.  However,  such  an  approach 
involves  solving  a  set  of  nonlinear  equations  which  are  sensitive  to  noise  in  the 
flow.  We  have  developed  a  new  algorithm  which  is  based  on  the  flow  field  and  its 
first  derivatives,  over  three  consecutive  frames,  under  the  assumption  that  the  3- 
D  motion  is  constant  over  time  and  object  surface  is  locally  planar.  This  algo¬ 
rithm  relies  on  solving  linear  systems  of  equations  and  has  proven  to  be  robust 
under  the  presence  of  noise. 

References: 

1.  K.  Wohn,  “Image  Deformation:  Its  Recovery  and  Utility,”  Proceedings  of 

International  Conference  on  Advances  in  Image  Processing  and  Pattern 

Recognition,  Pisa,  Italy,  December  1985. 


K.  Wohn  and  J.  Wu,  “Recovering  3-D  Structure  and  Motion  from  First- 
Order  Image  Deformation,”  manuscript. 


ID.3ProbabiIistic  Methods  in  Single  Frame  Picture  Processing.  R.W. 
Brockett,  J.  Clark,  and  A.  Mclvor,  Contracts  N00014-84-K-0465,  N00014- 
84-K-0504,  D AAG29-83-K-0027 ,  and  AFOSR-81-0054;  Research  Unit  9. 

Many  tracking  algorithms  approach  the  problem  in  two  stages.  The  first  is 
that  of  segmentation  and  the  second  is  that  of  nonlinear  filtering,  using  the  results 
of  segmentation  as  an  input.  Because  the  errors  made  in  the  segmentation  process 
impose  a  limit  on  the  performance  of  the  system  it  is  important  to  do  this  seg¬ 
mentation  as  accurately  as  possible.  In  his  thesis  [1]  Mclvor  has  put  the  problem 
of  edge  smoothing  in  a  filtering  context.  He  models  the  curve  edge  detector  pro¬ 
cess  with  arc  length  being  the  independent  variable.  This  process  is  thus  treated 
as  a  signal  and  the  output  of  a  crude  edge  detector  is  used  as  the  observation  pro¬ 
cess.  When  the  curve  is  taken  to  be  a  process  with  Gauss-Markov  curvature  the 
model  equations  are  bilinear.  This  situation  is  shown  to  give  good  resuls  in 
specific  cases. 

Reference 

1.  A.  Mclvor,  “Stochastically  Based  Vision  Algorithm,”  Ph.D  Thesis,  Harvard 
University,  June  1985. 
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PUBLICATIONS /PATENTS/PRESENTATIONS /HONORS 


a.  Papers  Submitted  to  Refereed  Journals  (and  not  yet  published) 

R.W.  Brockett,  “Robotic  Hand  With  Rheological  Surfaces,”  Proceedings  of 
the  1985  IEEE  Conference  on  Robotics  and  Automation,  St.  Louis,  MO. 


R.W.  Brockett  and  J.  Loncaric,  “The  Geometry  of  Compliance  Program¬ 
ming,”  MTNS  Proceedings  (to  appear). 

A.  Waxman  and  K.  Wohn,  “Contour  Evolution,  Neighborhood  Deformation 
and  Global  Image  Flow:  Planar  Surfaces  in  Motion,”  to  appear  in  Inti.  Jour¬ 
nal  of  Robotics  Research  4,  MIT  Press,  1985. 

K.  Wohn  and  J.  Wu,  “Recovery  of  3-D  Motion  of  Local  Planar  Surfaces  from 
Optical  Flow  Field,”  (in  progress). 

K.  Wohn,  “Image  Deformation:  Its  Recovery  and  Utility,”  to  appear  in  Inti. 
Conf.  Advances  in  Image  Processing  and  Pattern  Recognition,  Pisa,  Italy, 
December  1985. 

i.  Honors/Awards/Prizes 

A.M.  Mclvor,  “Stochastically  Based  Vision  Algorithms,”  Thesis,  Harvard 
University,  Division  of  Applied  Sciences,  1985. 

T.B.  Sanger,  “Stereo  Disparity  Computation  Using  Gabor  Filters,”  Thesis, 
Harvard  University,  Division  of  Applied  Sceinces,  1985. 

J.  Loncaric,  “Geometrical  Analysis  of  Compliant  Mechanisms  in  Robotics,” 
Thesis,  Harvard  University,  Division  of  Applied  Sciences,  1985. 


III.4Discrete  Event  Dynamic  Systems  Study.  Y.C.  Ho,  R.Suri,  X.  Cao,  J. 
Dille,  and  M.  Zazanis,  Contracts  N00014-84-K-0465,  N00014-79-C-0776, 
NSF-ECS-82- 13680,  and  NSF-CDR-85-001-08;  Research  Unit  10. 


The  overall  aim  of  the  researches  in  this  unit  is  to  develop  analysis,  optimiza 
tion,  and  control  techniques  for  man-rhade  systems.  Our  current  emphasis  is  on  a 
class  of  problems  called  discrete  event  dynamic  systems.  Examples  are 
computer/communication  networks,  automated  manufacturing  systems,  database 
under  concurrent  use.  to  name  just  three.  Our  earlier  work  in  the  previous  three 
year  funding  cycle  showed  that  a  new  technique  called  perturbation  analysis  has 
considerable  promise.  The  goal  here  is  to  attempt  the  extension  and  generaliza¬ 
tion  of  this  technique  to  a  broad  class  of  discrete  event  systems  in  analogy  to  the 
past  successes  we  had  with  the  more  familiar  dynamic  systems  governed  by 
differential  equations. 


Perturbation  analysis  provides  a  very  efficient  way  to  estimate  the  sensitivi¬ 
ties  of  system  performances  with  respect  to  parameters  for  discrete  event  systems. 
Algorithms  have  been  developed  for  various  queueing  networks  to  obtain  these 
estimates  based  on  only  a  single  sample  path.  Extensive  experiments  show  that 
these  estimates  are  quite  accurate.  Recently,  some  significant  accomplishment 
concerning  the  convergence  properties  of  these  estimates  has  been  achieved.  The 
convergence  theory  shows  that  the  estimates  obtained  using  perturbation  analysis 
are  strongly  consistent,  hence  justifying  the  previous  algorithms. 

Two  main  performance  measures  of  a  closed  Jackson  network,  the  system 
throughput  and  the  mean  sojourn  times  of  a  customer  staying  in  a  particular 
server  or  the  system,  are  studied.  It  is  proved  that  perturbation  analysis  provides 
a  strongly  consistent  estimate  of  the  derivative  of  the  expected  throughput  with 
respect  to  the  mean  service  time  in  a  finite  period.  The  elasticity  of  the 
throughput  with  respect  to  the  mean  service  time  converges  in  the  mean  to  the 
elasticity  of  the  throughput  in  steady  state  as  the  number  of  customers  served 
goes  to  infinity.  Similar  results  hold  for  the  sensitivity  of  the  mean  sojourn  times. 
We  also  developed  a  set  of  new  equations  based  on  perturbation  analysis  which 
calculates  theoretically  the  derivatives  of  throughput  with  respect  to  the  mean 
service  time.  Numerical  solutions  show  that  the  results  are  exactly  the  same  as 
those  obtained  by  the  well-known  Jackson  formulas. 

The  sensitivity  of  mean  sysem  times  of  a  customer  to  a  parameter  of  the 
arrival  or  service  distribution  in  a  M/G/l  queue  is  also  studied.  It  is  shown  that, 
perturbation  analysis  gives  asymptotically  unbiased  and  strongly  consistent  esti¬ 
mates  of  this  sensitivity. 

Meanwhile,  the  mean  square  error  of  the  perturbation  analysis  estimates  is 
examined.  The  error  is  compared  with  that  of  estimates  obtained  by  conventional 
methods.  The  results  establish  the  asymptotic  superiority  of  perturbation 


analysis  over  conventional  methods. 

The  convergence  property  of  perturbation  analysis  estimates  for  more  general 
systems  is  now  under  investigation. 


ANNUAL  REPORT  OF 

PUBLICATIONS/P  ATENTS/PRESENT/.TIONS/HONORS 


a.  Papers  submittted  to  Refereed  Journals  (and  not  yet  published) 

Y.C.  Ho  and  W.B.  Gong,  “A  Note  on  Filtering  in  Queueing  Networks  and 
Discrete  Event  Dynamic  Systems,”  submitted  to  IEEE  Trans,  on  Automatic 
Control. 

X.  Cao,  “First  Order  Perturbation  Analysis  of  Multiclass  Queueing  Net¬ 
works,”  to  be  presented  at  the  1985  IEEE  Conference  on  Decision  and  Con¬ 
trol,  submitted  to  Performance  Evaluation. 

Y.  Dallery  and  R.  Suri,  “Approximate  Disaggregation  and  Performance 
Bounds  for  Queueing  Networks  wih  Multiple-Server  Stations,”  submitted  to 
PERFORMANCE,  1986. 

X.  Cao,  “On  the  Sample  Functions  of  Queueing  Networks  with  Application 
to  Perturbation  Analysis,”  submitted  to  Operations  Research. 

X.  Cao  and  Y.  Dallery,  “An  Operational  Approach  to  Perurbation  Analysis 
of  Closed  Queueing  Networks,”  Manuscript,  Harvard  University,  September 
1985,  submitted  to  1986  ACC  Conference. 

M.  Zazanis  and  R.  Suri,  “Comparison  of  Perturbation  Analysis  with  Conven¬ 
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IV.  ELECTROMAGNETIC  PHENOMENA 


Personnel 

Prof.  T.T.  Wu  Mr.  M.F.  Brown 

Prof.  R.W.P.  King  Ms.  M.  Owens 

Dr.  J.T.  DeBettencourt  Mr.  W.-Y.  Pan  (until  7/31/85) 

Dr.  J.M.  Myers  (after  11/1/85)  Ms.  B.H.  Sandler 

Research  in  the  area  of  electromagnetic  radiation  is  directed  toward  the  solu¬ 
tion  of  practical  problems  through  the  complete  understanding  of  the  underlying 
physical  phenomena.  This  involves  the  coordinated  application  of  modern  analyti¬ 
cal,  numerical,  and  experimental  techniques  and  the  use  of  high-speed  computers 
and  precision  instrumentation.  Application  is  also  made  of  modeling  techniques 
and  the  principle  of  similitude.  Most  practically  significant  problems  in  the  area 
are  sufficienly  complicated  that  extensive  computation  and  measurement  are  often 
required  to  justify  approximations  that  are  usually  necessary.  Where  possible, 
general  formulas  are  obtained  and  verified  experimentally  so  that  the  phenomenon 
under  study  can  be  understood  physically  in  analytical  form  and  not  just  as  a  set 
of  numbers. 

The  researches  jare  concerned  primarily  with  the  properties  of  antennas  and 
arrays  and  of  the  electromagnetic  fields  they  generate  in  various  practically 
important  environments  that  lead  to  difficult  problems  with  complicated  boun¬ 
dary  conditions.  Examples  include  dipoles,  traveling-wave  antennas  and  arrays, 
crossed  dipoles,  and  loops  near  the  boundary  between  two  media  such  as  air  and 
the  earth  or  sea,  or  the  oceanic  crust  and  sea  water;  the  scattering  of  electromag¬ 
netic  waves  from  buried  or  submerged  objects;  the  properties  of  lateral  elec¬ 
tromagnetic  waves  and  pulses;  the  field  inside  and  currents  and  charges  induced 
(  hstacles  within  an  EMP  simulator;  arrays  of  antennas  along  curved  lines;  and 


solitary  electromagnetic  pulses  with  slow  rates  of  decay. 


IV.lNumerical  and  Analytical  Determination  of  the  Fields  of  Antennas 
near  an  Interface  Between  Two  Half-Spaces  with  Significantly 
Different  Wave  Numbers.  T.T.  Wu,  R.W.P.  King,  B.H.  Sandler,  and  M. 
Owens,  Contract  N00014-84-K-0465  and  N00014-79-C-0419;  Research  Unit 
11. 

The  recently  derived  sets  of  simple  and  accurate,  analytical  formulas  for  the 
complete  electromagnetic  fields  of  horizontal  [1-3]  and  vertical  [4]  electric  dipoles 
on  or  near  the  boundary  between  two  half-spaces  have  provided  new  insight  into 
lateral  propagation  in  general,  and  more  flexible  tools  for  a  variety  of  specific 
applications.  Examples  of  their  usefulness  are  presented  here  and  under  topic  #2. 

An- examination  of  the  new  formulas  shows  that  surface  waves  have  unusual 
properties  that  make  them  valuable  tools  for  geophysical  prospecting  [5,6],  They 
show  that  over  all  distances  and  at  all  frequencies  a  signal  reaching  a^ receiving 
dipole  in  the  salt  water  near  the  sea  floor  has,  in  effect,  proceeded  downward  from 
the  transmiting  dipole  in  the  sea  water  into  the  rock,  then  radially  outward  in  the 
rock,  and  finally  upward  and  back  into  the  sea  water  to  the  receiving  dipole.  The 


decrease  in  amplitude  and  change  in  phase  with  radial  distance  depend  primarily 
on  the  electrical  properties  of  the  rock.  With  a  suitable  choice  of  frequencies  and 
radial  distances  p  between  transmitting  and  receiving  dipoles,  the  received  signal 
can  have  the  simple  form  exp(-ap)exp(i0p)/p  where  a  is  the  attenuation  constant 
and  P  is  the  real  wave  number  of  the  rock.  These  important  facts  are  the  basis 


for  the  development  of  methods  for  the  determination  of  the  conductivity  and  per¬ 
mittivity  of  the  rock  forming  the  sea  floor  by  means  of  measurements  made  in  the 
adjacent  sea  water.  Since  lateral  waves  do  not  penetrate  deeply  into  the  ocean 


floor,  they  are  useful  primarily  in  determining  the  average  conductivity  and  per- 
mitivity  of  the  material  fairly  close  to  the  boundary.  For  this  a  two-layer  theory 
may  be  adequate. 
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The  new  formulas  have  also  been  used  [5,6]  to  assist  in  the  interpretation  of 
available  measurements  of  the  conductivity  of  the  sea  floor.  It  appears  that  the 
measured  data  would  be  well  satisfied  by  a  half-space  model  of  the  lithosphere 
that  is  homogeneous  but  one-dimensionally  anisotropic  in  the  conductivity.  Such 
a  model  involves  no  reflections  from  deeper,  widely  spaced  horizontal  layers.  An 
analytical  study  by  W.-Y.  Pan  of  the  one-dimensionally  anisotropic  half-space  has 
been  completed  which  agrees  with  the  available  measured  data.  It  is  discussed 
under  topic  No.  3. 

Currently  the  new  formulas  are  being  used  to  determine  certain  properties  of 
lateral  waves.  As  a  first  step  in  the  process,  explicit  formulas  have  been  obtained 
[7,8]  for  the  complete  field  in  Region  2  (e.g.,  air,  lithospheric  rock)  of  a  horizontal 
electric  dipole  located  near  the  boundary  in  Region  1  (e.g.,  -  .It  water).  Expres¬ 
sions  for  the  field  in  Region  2  are  required  in  order  to  determine  such  properties  as 
the  locus  of  the  Poynting  vector,  the  maximum  depth  of  penetration  of  the  lateral 
wave  into  Region  2,  and  the  fraction  of  power  associated  with  the  lateral  wave. 
The  study  is  continuing  with  the  determination  of  the  field  in  Region  2  of  the 
vertical  electric  dipole. 
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IV.2Theoretical  Study  of  Scattering  from  Buried  or  Submerged 

Objects.  R.W.P.  King,  T.T.  Wu  and  W.-Y.  Pan,  Contract  N00014-84-K- 
0465;  Research  Unit  11. 

The  detection  and  localization  of  buried  or  submerged  objects  or  inhomo¬ 
geneities  in  a  dissipative  half-space  from  measurements  of  a  scattered  electromag¬ 
netic  field  are  of  current  interest.  Earlier  investigations  have  been  concerned  pri¬ 
marily  with  plane  waves  incident  on  the  surface  of  the  earth  from  the  air  above  at 
an  arbitrary  angle.  The  present  study  deals  with  the  incident  fields  generated  by 
horizontal  and  vertical  antennas  located  on  the  surface  of  the  earth.  The  first 
step  in  the  analysis  is  the  determination  of  this  incident  field.  The  complete  field 
of  such  antennas  can  now  be  expressed  in  simple  and  accurate  form  with  the  avai¬ 
lability  of  new  formulas  (see  topic  No.  1,  refs.  1-4)  that  are  excellent  approxima¬ 
tions  of  the  rigorous  general  integrals.  In  the  first  of  two  papers  [l],  the  com¬ 
ponents  of  the  lateral-wave  field  at  points  near  the  surface  in  air  and  at  all  points 
in  the  earth  have  been  determined  when  the  source  is  a  horizontal-wire  antenna  of 
finite  length  located  a  small  height  over  the  earth  or  salt  water.  Both  antennas 
with  open  ends  and  antennas  with  ends  grounded  by  vertical  conductors  are 
treated. 

Once  the  incident  field  is  known,  the  scattered  field  can  be  investigated.  The 
first  scattering  object  studied  is  the  horizontal  insulated  conductor  [2].  since  this 
provides  a  scattered  field  that  is  readily  calculated  and  combined  with  the 
incident  field  to  obtain  a  total  field  that  exhibits  interference  phenomena.  The 


incident  lateral-wave  field  is  maintained  by  a  horizontal  dipole  in  air  just  above 
the  surface  of  the  earth,  as  described  above.  The  total  field  (incident  plus  scat¬ 
tered  fields)  is  to  be  observed  also  in  the  air  just  above  the  surface.  The  current 
induced  in  the  insulated  wire  by  the  lateral-wave  field  incident  from  the 
transmitter  is  evaluated  first  and  then  the  reradiated  field  just  above  the  surface 
is  calculated.  Two  orientations  are  considered:  The  field  at  the  center  of  the 
buried  conductor  is  (1)  parallel  to  Elf(p,0,z),  and  (2)  parallel  to  El4,{p,n/2,z).  The 
incident,  scattered  and  total  fields  in  an  area  above  the  conductor  have  been 
evaluated  for  the  conditions  of  a  laboratory  model.  It  is  clearly  demonstrated 
that  a  significant  change  in  the  electric  field  over  the  volume  occupied  by  the 
scattering  object  is  observed,  and  that  this  is  sufficiently  localized  to  permit  an 
accurate  bounding  of  the  area  above  the  object. 

Only  a  few  scattering  problems  involving  very  simply  shaped  targets  can  be 
solved  analytically.  In  a  separate  study  by  W.-Y.  Pan,  analytical  expressions 
have  been  obtained  [3]  for  the  scattered  field  from  a  buried,  perfectly  conducting 
rectangular  plate  placed  parallel  to  the  interface  between  air  and  earth.  First,  the 
current  is. expressed  in  terms  of  several  Tchebichef  polynomials  of  the  second  kind. 
The  scattered  field  is  the  sum  of  the  fields  that  are  generated  separately  by 
equivalent  horizontal  dipoles  and  multipoles.  The  moments  of  the  dipoles  and 
multipoles  are  obtained  from  the  coefficients  of  the  current  expansion.  These 
coefficients  have  been  calculated  over  a  frequency  range  from  13  to  160  kHz  and  a 
range  of  plate  sizes  from  2  to  20,000  square  meters  in  area. 
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3.  W.-Y.  Pan,  “The  Scattered  Field  from  a  Buried,  Perfectly  Conducting  Plate 

in  a  Lossy  Earth,”  Electromagnetics,  submitted  for  publication. 


IV.3Theoretical  Study  of  the  One-Dimensionally  Anisotropic  Half- 

Space.  W.-Y.  Pan,  Contract  N00014-84-K-0465;  Research  Unit  11. 

Anisotropy  in  conductivity  is  found  in  various  stratified  media  including 
alternating  layers  of  dense  rock  with  low  conductivity  and  less  dense  rock  with 
higher  conductivity.  In  such  media  the  conductivity  transverse  to  the  bedding 
surfaces  is  always  smaller  than  that  along  the  bedding  surfaces.  A  model  of  hor¬ 
izontally  stratified  rock  in  the  form  of  relatively  thin  layers  that  alternate  with 
higher  and  lower  conductivities  is  well  approximated  by  a  half-space  that  is  homo¬ 
geneous  but  anisotropic  in  conductivity,  with  the  conductivity  <xz  =  <rL  perpendic¬ 
ular  to  the  bedding  surfaces  different  from  the  conductivity  <r1  =  av  =  <rT  parallel 
to  these  surfaces.  Both  aL  and  <rT  are  constant.  An  analytical  study  of  the 
one-dimensionally  anisotropic  half-space,  beginning  with  Maxwell’s  equations  and 
paralleling  the  analysis  for  the  isotropic  half-spaces  (see  topic  No.  1,  refs.  1-3),  has 
been  completed  [l]  and  compared  with  recent  measurements  of  the  oceanic  crust 
[2].  It  confirms  the  validity  of  the  half-space  model  with  o,  =  oy  =  0.004S/m  and 
<t:  =  0.002S/m,  and  explains  the  different  behavior  observed  for  the  E and  Ey 
components.  Reflections  from  deeper,  widely  spaced  horizontal  layers  are  not 
required  to  explain  the  observed  data.  When  <rL  =  <rT,  the  new  formulas  coincide 
with  those  of  Wu  and  King  for  an  isotropic  medium.  When  <rL  *  <rT,  the  lateral 
wave  part  of  the  field  is  separated  into  two  terms  of  which  the  first  is  of  electric 
type  with  the  wave  number  k2  =  kL,  and  the  second  is  of  magnetic  type  with 


A?2  38  kT. 
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IV.4Antennas  in  Material  Media  near  Boundaries  and  Their  Applica¬ 
tion.  R.W.P.  King,  contract  N00014-84-K-0465;  Research  Unit  11. 

Antennas  are  located  on  or  near  the  boundary  between  two  electrically 
different  media  for  many  purposes  including  communication,  remote  sensing,  and 
geophysical  exploration.  Their  properties  and  the  field  they  generate  are  strongly 
influenced  by  the  proximity  of  the  boundary.  Some  of  the  effects  of  a  boundary 
are  conveniently  represented  in  terms  of  a  suitable  image.  Others,  like  the 
lateral-wave  fields,  are  in  addition  to  image  effects.  This  is  well  understood  for 
vertical  antennas  erected  over  the  earth  or  sea  [l],  and  for  horizontal  wave  anten¬ 
nas  on  the  surface  of  the  earth  [2].  In  a  new  study  [3,4]  the  properties  of  antennas 
embedded  in  a  material  medium ’near  a  boundary  have  been  determined.  In  the 
first  paper  [3],  the  theory  of  the  bare  metal  dipole  antenna  has  been  developed  in 
terms  of  two  specific  applications:  the  subsurface  antenna  for  communication 
between  two  points  below  the  surface  of  the  earth,  and  the  seafloor  antenna  for 
use  in  studies  of  the  conductivity  of  the  oceanic  crust.  In  the  first  case,  the 
antenna  is  embedded  in  a  medium  like  soil,  sand,  ice,  fresh  or  salt  water  near  its 
boundary  with  air;  a  frequency  range  from  30  to  300  kHz  is  useful.  In  the  second, 
the  antenna  is  located  in  the  salt  water  near  the  oceanic  crust  under  the  sea  floor; 
in  this  case  very  low  frequencies  in  the  range  from  0.3  to  3  Hz  are  needed  to 
penetrate  sufficiently  into  the  lithosphere.  In  both  locations,  the  transmitting  and 
receiving  antennas  are  located  in  a  region  of  large  wave  number  (Region  1) 


adjacent  to  a  region  with  a  much  smaller  wave  number  (Region  2). 

The  principal  differences  between  the  antennas  near  the  air  surface  and  those 
near  the  sea  floor  are  a  consequence  of  the  large  difference  in  frequency.  At  the 
very  low  frequencies  used  on  the  sea  floor,  the  internal  impedance  of  the  copper 
conductor  is  important,  whereas  it  is  negligible  at  the  higher  frequencies  used  near 
the  air  surface.  This  means  that  whereas  the  bare  metal  antenna  is  highly 
efficient  in  transferring  power  from  the  antenna  to  the  exterior  media  (Regions  1 
and  2)  at  the  higher  frequencies,  it  is  subject  to  an  unavoidable  large  internal 
power  loss  in  heating  the  antenna  at  the  low  frequencies.  The  properties  of  bare 
copper  antennas  in  the  two  applications  are  summarized  [3].  These  include  the 
distribution  of  current,  the  driving-point  admittance  and  impedance,  and  the 
effective  length. 

When  the  ambient  medium  in  which  a  bare  metal  antenna  is  embedded  is 
conducting,  radial  currents  from  the  antenna  into  the  medium  decrease  the  axial 
current.  This  means  that  often  only  electrically  relatively  short  antennas  are  use¬ 
ful  with  correspondingly  small  effective  lengths.  Radial  currents  can  obviously  be 
eliminated  by  enclosing  the  metal  cylinder  in  a  dielectric  sheath.  This  greatly 
changes  the  properties  of  the  antenna.  In  the  second  paper  [4],  the  properties  of 
the  insulated  antenna  in  a  dissipative  medium  near  a  boundary  are  reviewed.  The 
important  effects  on  antennas  of  being  embedded  in  media  with  different  proper¬ 
ties  over  a  wide  range  of  frequencies  are  treated  spedifically  in  terms  of  three 
examples:  (l)  A  traveling-wave  antenna  in  the  earth  below  air  at  /  ~  0.1  MHz  for 
subsurface  communication;  (2)  an  antiresonant  antenna  on  the  sea  floor  at  /  ~  l 
kHz;  (3)  an  electrically  short  antenna  on  the  sea  floor  at  /  ~  l  Hz.  The  last  two 
are  both  for  geophysical  measurements.  In  each  case  the  current  distribution, 
impedance,  input  power,. effective  length,  and  electric  field  in  the  range  of  interest 
are  determined. 


It  is  shown  that  the  insulated  conductor  with  its  various  possible  termina¬ 
tions  (including  especially  bare  monopoles)  is  a  highly  effective  antenna  when 
embedded  in  a  general  dielectric  or  conducting  region  with  a  wave  number  that  is 
much  larger  in  magnitude  that  the  wave  number  of  the  insulating  sleeve.  Depend¬ 
ing  on  the  ratio  6/a  of  the  radius  of  the  insulating  sleeve  to  that  of  the  conduc¬ 
tor,  the  insulated  conductor  acts  more  like  an  antenna  for  transferring  power  to 
the  ambient  medium  (6/a  ~  l)  or  a  transmission  line  for  supplying  power  to  the 
terminations  (6/a  >  10).  At  frequencies  of  1  kHz  and  higher,  the  power  loss  in  the 
internal  resistance  of  the  conductor  is  negligible  and  all  power  is  either  radiated  or 
transferred  to  the  terminations  without  loss.  At  frequence  ^  below  100  Hz,  the 
power  loss  in  the  internal  resistance  becomes  dominant  aid  i  he  radiated  power 
transferred  directly  to  the  ambient  medium  negligible  small.  In  this  case,  the 
insulated  conductor  is  a  highly  inefficient  transmission  line  feeding  the  termina¬ 
tions. 
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IV.5Experimental  Study  of  Lateral-Wave  Propagation  in  Layered 

Media.  W.-Y.  Pan,  Contract  N00014-84-K-0465;  Research  Unit  11. 

Current  interest  in  the  possibility  of  investigating  the  electrical  properties  of 
the  earth  beneath  the  ocean  by  electromagnetic  means  naturally  leads  to  the 
problem  of  wave  propagation  in  layered  media.  A  realistic  model  of  the  earth’s 
crust  typically  includes  a  layer  of  sediment  followed  by  many  layers  of  rock  of 
various  types.  Lateral-wave  propagation  for  such  cases  is  much  more  complicated 
since  a  different  wave  can  travel  in  each  of  the  rock  regions  and,  at  the  same  time, 
each  of  the  waves  is  influenced  by  the  neighboring  regions.  The  recent  theoretical 
work  of  Dunn  [l,2]  on  a  three-layer  model  provides  generalized  formulas  for 
lateral  waves  generated  by  a  horizontal  electric  dipole  with  unit  moment  in  sea 
water  (Region  1)  near  its  boundary  with  a  medium  (Region  2)  of  arbitrary  thick¬ 
ness  and  having  electrical  properties  intermediate  to  those  of  the  sea  water  and 
the  underlying  rock  (Region  3).  His  three-layer  solution  looks  very  much  like  the 
Wu  and  King  two-layer  solution.  The  differences  are:  (1)  there  is  a  sediment  term 
which  represents  the  contribution  from  the  wave  traveling  in  the  sediment  layer; 
at  distances  and  frequencies  of  interest  it  is  negligible;  (2)  there  is  a  constant  A 
which  can  change  the  magnitudes  of  the  fields;  and  (3)  the  Fresnel  integral  now 
depends  on  a  complicated  constant  u. 

During  this  reporting  period,  experiments  corresponding  to  this  three-layer 
configuration  have  been  performed  in  the  laboratory  using  a  model  tank  filled  with 
distilled  water.  The  apparatus  also  included  a  generator,  vector  voltmeter,  . 
transmitting  antenna,  and  automated  movable  receiving  antenna.  The  octagonal 
tank  measures  8.5  feet  in  diameter  and  2.5  feet  in  height.  The  transmitting 
antenna  is  located  at  the  center  of  the  tank;  the  receiving  antenna  can  move  radi¬ 
ally  outward  from  the  transmitter  on  a  track.  Both  of  the  antennas  are  horizon¬ 
tal  insulated  linear  dipoles,  oriented  in  the  radial  direction.  In  order  to  have  the 
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radial  range  extend  a  few  wavelengths,  the  frequency  was  chosen  to  be  1.5  GHz. 

On  the  water  surface  there  is  a  slab  supported  by  four  plastic  columns.  This 
slab  is  also  octagonal  and  a  little  smaller  then  the  tank  in  size.  It  consists  of  one 
or  more  sheets  of  homosote  adhered  to  a  slab  of  styrofoam,  16  cm  thick.  The 
layers  are  held  together  by  plastic  bolts  and  nuts  to  assure  flatness.  Initially  the 
composite  slab  was  wrapped  with  very  thin  plastic  film  to  keep  the  homosote  dry. 
After  the  field  had  been  measured  with  the  dry  homosote,  the  plastic  film  was 
removed  and  the  homosote  was  lowered  into  the  water  for  three  days  to  absorb 
water.  In  the  process,  the  thickness  of  each  sheet  of  homosote  increased  from  1.35 
cm  to  1.6  cm.  By  weighing  a  small  sample  of  the  homosote,  it  was  estimated  that 
the  completely  wet  homosote  slab  included  18  percent  water.  Since  the  conduc¬ 
tivity  and  relative  dielectric  constant  of  the  styrofoam  are  nearly  the  same  as 
those  of  air,  it  is  unnecessary  to  distinguish  the  styrofoam  from  the  air.  The 
electrical  properties  of  the  homosote  were  estimated  to  be:  =  2.7  and 

c ?■ 2dry  =  0.01352  S/m  for  the  dry  homosote;  <£*‘  =  16.61  and  <x2wel  =  0.0993  S/m  for  t  he 
wet  homosote.  When  the  homosote  is  wet,  the  electrical  properties  in  the  three 
regions  approximately  satisfy  the  condition  l/t,l  >  31/bJ  >  3|jtjl,  as  required  in  Dunn's 
theory. 

The  radial  electric-field  component  E lp  was  measured  with  one,  two,  and 
four  sheets  of  homosote  in  both  the  dry  and  wet  states  and  a  paper  has  been  writ¬ 
ten  comparing  the  measured  results  with  theory  [3].  The  following  conclusions 
have  been  drawn:  (1)  The  measurements  for  the  wet  slab  agree  quite  well  with 
the  results  calculated  using  Dunn’s  theory.  The  field  at  the  receiving  antenna  is  a 
lateral  wave  which  is  generated  by  the  transmitting  antenna,  travels  up  through 
the  middle  layer  of  homosote,  then  radially  forward  in  the  air,  and  finally  down 
through  the  middle  layer  into  the  water  and  to  the  receiving  antenna.  (2)  In  the 
near  range,  the  field  with  the  wet  homosote  slab  will  be  larger  than  the 


corresponding  value  without  the  middle  layer.  (3)  If  the  thickness  of  the  middle 
layer  is  greater  then  one  wavelength  in  the  medium,  the  amplitude  of  the  wave 
propagating  in  the  middle  layer  is  not  small  and  the  wave  will  interfere  with  the 
lateral  wave  propagating  in  the  air  in  a  range  of  distance.  (4)  When  the  electrical 
properties  of  the  middle  layer  are  not  very  different  from  those  of  air  (as  is  the 
case  for  dry  homosote),  the  field  is  similar  to  that  of  the  two-layer  problem. 
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IV.6LateraI-Wave  Propagation  and  Modeling  of  the  Lithosphere.  M.F. 

Brown,  R.W.P.  King,  T.T.  Wu,  and  J.T.  deBettencourt,  Contracts  N00014- 

79-C-0419  and  N00014-84-K-0465;  Research  Unit  11. 

Completion  of  a  detailed  charting  of  the  transient  response  of  vertical  linear 
monopoles  above  a  metal  ground  plane  has  been  attained.  Groups  of  transient 
waveforms  have  been  plotted  as  a  function  of  antenna  length  (4  cm,  13.5  cm,  and 
from  20  to  140  cm  in  various  increments)  and  source-receiver  interantenna  separa- 
tions  (from  20  cm  to  180  cm  in  20-cm  increments).  A  quite  clean  received  pulse 
(one  which  is  very  similar  in  form  to  Gaussian  output  of  pulse  generator)  has  been 
measured  using  a  140-cm  tall  monopole,  and  chosen  as  a  reference  pulse.  A  linear 
telescoping  “TV”  antenna,  when  extended  beyond  140  cm  above  its  underlying 
metal  ground  plane,  has  been  shown  to  produce  little  additional  improvement  in 
received  waveform  relative  to  the  reference  pulse. 


The  aforementioned  metal  ground  plane,  supported  by  the  walls  of  a  shallow 
(~  1-J./2  ft)  wading  pool,  has  been  removed  to  expose  an  underlying  region  of  salt 
water  with  ocean  constitutive  parameters  (<r  =  3.5S/m,  tT  =  80).  Both  short  (4.0 
cm)  and  tall  (140  cm)  vertical  monopole  antennas  have  been  placed  in  a 
configuration  chosen  to  replicate  the  ground  plane  studies,  with  one  exception.  A 
small  metal  circular  disk  (~  4  cm  in  diameter)  is  used  as  a  conducting  plate  at  the 
base  of  the  antenna  where  it  contacts  the  salt  water.  The  conducting  plate  is  sup¬ 
ported  bv  vertical  rods  which  extend  to  the  floor  of  the  wading  pool.  The  rods 
and  their  base  at  the  wading  pool’s  floor  produce  a  negligible  effect  upon  the  field 
measured  in  air  by  the  vertical  monopole  receiver  whose  bottom  tip  is  ~  3  cm 
above  the  salt-water  surface. 

The  measured  waveforms  at  various  source-receiver  separations  (from  20  cm 
to  140  cm  in  20-cm  increments)  for  vertical  monopoles  above  salt  water  have  been 
shown  to  be  quite  similar  to  those  for  dipoles  above  the  metal  ground  plane  of  the 
preceding  series  of  data.  The  center  frequency  of  the  received  pulses  has  been 
measured  to  be  ~  500  MHz.  The  pulses  undergo  a  decrement  in  peak  amplitude 
and  a  small  degree  of  spreading,  with  radial  distance  from  the  source  monopole. 
The  Fourier  transforms  of  all  waveforms  have  been  concurrently  obtained  with 
plots  of  E{t)  vs.  t,  and  graphed  using  the  waveform  digitizer  (Tektronix 
WP2110).  The  Fourier  transforms  of  the  pulses  which  propagate  over  salt  water 
indicate  that  the  high-frequency  components  of  the  received  waveform  decrease 
with  source-receiver  separation.  A  paper  on  the  transient  response  of  the  vertical 
monopole  over  a  metal  ground  plane  and  over  salt  water  has  been  presented  [  1  ] . 

Studies  are  currently  underway  to  categorize  the  waveform  due  to  submerged 
horizontal  pulsed  sources  in  the  form  of  an  insulated  antenna  and  broadside  arrays 
of  such  elements.  Although  the  theory  of  insulated  antennas  for  use  in  CW  work 
indicated  that  they  are  quite  frequency-dependent,  a  choice  of  design  associated 


with  the  center  frequency  of  the  input  pulse  has  yielded  a  very  strong  received 
waveform  (using,  as  in  all  measurements  referred  to  above,  a  vertical  monopole 
receiver  in  a  noise- reducing  configuration).  The  received  pulse  is  in  the  familiar 
form  of  a  Gaussian  field  E(t)  vs  t,  and  the  initial  2-ns  wide  input  pulse  undergoes 
spreading  and  decrease  in  amplitude  similar  to  that  observed  for  vertical  mono¬ 
poles  above  salt  water. 

Reference: 

1.  M.F.  Brown,  “The  Pulsed  Lateral-Wave  Field  of  a  Vertical  Monopole  in  Air,’" 
presented  at  the  1985  USNC/URSI  Commission  F  Meeting,  held  at  U.  Mass., 
Amherst,  Oct.  1985. 


IV.7Fields  and  Currents  and  Charges  on  Obstacles  in  a  Parallel-Plate 
Simulator  at  Selected  Frequencies  and  with  Pulse  Excitation.  W.-Y. 
Pan,  R.W.P.  King,  and  T.T.  Wu,  Contract  N00014-84-K-0465;  Research  Unit 
11. 

An  important  first  step  in  understanding  the  electromagnetic  field  which 
penetrates  into  the  interior  of  a  rocket  or  aircraft  through  apertures  in  the  metal 
skin,  when  this  is  exposed  to  an  electromagnetic  pulse  from  a  nuclear  explosion,  is 
the  determination  of  the  distributions  of  current  and  charge  on  the  outside  of  the 
metal  surface.  For  simplicity  the  shape  of  the  rocket  can  be  approximated  by  a 
finite  tubular  cylinder.  The  distributions  of  current  and  charge  induced  on  such  a 
cylinder  by  a  continuous-wave  field  have  been  studied  both  theoretically  [l]  and 
experimentally  in  the  Harvard  parallel-plate  simulator  [2].  A  corresponding  exper¬ 
imental  investigation  of  the  distributions  of  current  and  charge  induced  directly 
on  the  cylinder  by  an  electromagnetic  pulse  has  now  been  completed  [3].  Fourier 
transforms  of  the  observed  pulse  shapes  are  taken  in  order  to  obtain  distributions 
of  a  series  of  frequency  components  for  the  current  and  charge  on  the  surface. 
These  are  then  compared  with  theoretical  distributions  calculated  using  Kao’s 
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method  [4].  The  measured  data  are  found  to  be  in  basically  good  agreement  with 
those  theoretical  distributions  induced  by  the  sum  of  the  incident  and  reflected 
pulses.  Observed  differences  between  measured  and  theoretical  distributions  for 
certain  frequency  components  of  the  transverse  current  I,  are  probably  due  to 
coupling  between  the  cylinder  and  the  simulator  top  plate.  A  quantiative  analysis 
of  the  interaction  between  the  cylinder  and  the  simulator  is  difficult  and  requires 
further  study. 

The  measured  current  pulse  sequences  induced  on  the  vertical  cylinder 
described  above  have  also  been  compared,  in  a  separate  paper  [5],  with  theoreti¬ 
cally  determined  pulse  shapes  of  current  induced  on  an  infinitely  long  conducting 
cylinder  by  an  E-polarized  electromagnetic  pulse.  For  an  ^-polarized  incident 
electromagnetic  pulse  with  a  nonzero  direct  component,  the  axially  induced 
current  pulse  has  an  infinite  direct  component  and  the  Fourier  integration  is 
divergent.  In  order  to  avoid  this  difficulty,  the  Laplace  transform  for  the  variable 
t  is  taken  instead  of  the  Fourier  transform.  After  the  inverse  Laplace  transform 
has  been  taken,  the  pulse  shapes  of  the  axial  induced  current  can  be  expressed  as 
the  sum  of  only  a  few  terms.  The  measured  pulse  sequences  for  the  finite  cylinder 
closely  resemble,  in  its  main  features,  those  of  the  infinite  cylinder  but  with  a 
series  of  small  reflected  pulses  added. 

The  experimental  and  theoretical  investigation  of  the  rhombic  simulator 
under  pulse  excitation  [6,7]  indicated  that  the  electromagnetic  pulse  in  the  work¬ 
ing  volume  of  the  simulator  is  excited  by  the  current  in  the  wires  near  the  source 
rather  than  that  in  the  wires  near  the  center  of  the  simulator.  This  observation  is 
somewhat  different  from  that  observed  in  the  frequency  domain.  The  same  ques¬ 
tion  must  be  asked  for  the  parallel-plate  simulator:  Does  the  triangular  plate  near 
the  source  also  play  the  most  important  role  in  forming  the  electromagnetic  pulse 
in  the  working  volume  of  the  parallel-plate  simulator?  To  answer  this  question. 
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the  electromagnetic  pulse  has  been  evaluated  in  terms  of  the  currents  on  the 
plates  and  the  roles  of  currents  on  different  parts  of  the  plates  have  been 
analyzed.  A  paper  has  been  written  which  clarifies  the  excitation  mechanism  of 
the  parallel-plate  EMP  simulator  [8].  It  is  found  that  if  the  length  of  the  parallel 
plate  is  short  (i.e.,  close  to  the  height  h),  the  currents  on  the  parallel  plate  contri¬ 
bute  almost  nothing  to  the  electromagnetic  pulse. 
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IV.8Advanced  Introduction  to  Electromagnetic  Theory  and  Its  Applica¬ 
tion.  R.W.P.  King  and  M.  Owens,  and  S.  Prasad  Hinchey  of  Northeastern 
University,  Contract  N00014-K-0465;  Research  Unit  11. 


A  new  introductory  book  on  electromagnetic  theory  and  its  applications  has 
been  published  by  Prentice-Hall  [l].  The  foundations  of  electromagnetism  are 
presented  in  the  first  six  chapters  in  a  manner  following  that  in  King’s  Fundamen¬ 
tal  Electromagnetic  Theory-- which  is  now  out  of  print.  Applications  begin  with  a 
chapter  on  the  scattering  and  diffraction  of  plane  waves  by  a  half-plane  and  con¬ 
tinue  with  chapters  on  antennas,  electric  circuit  theory,  transmission-line  theory, 
the  insulated  antenna,  the  theory  of  metal  waveguides,  waves  along  dielectric  rods 
with  reference  to  optical  fiber  transmission,  and  electromagnetic  surface  waves 
along  boundaries.  The  presentation  stresses  the  physical  interpretation  of  the 
mathematical  symbolism  and  the  basic  unity  of  the  theory.  The  applications  are 
arranged  with  increasing  difficulty  as  logical  continuations  of  the  theory.  They 
deal  specifically  with  problems  of  continuing  and  current  interest. 
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IV.OSolitary  Electromagnetic  Pulses  with  a  Slow  Rate  of  Decay.  T.T. 
Wu,  R.W.P.  King,  and  J.M.  Myers,  Contracts  N00014-84-K-0465  and 
N00014-85-K-0616  and  DOE  Grant  DE-FG02-84ER40158;  Supplementary 
Research,  Unit  11. 


In  a  recent  paper  by  J.M  Brittingham  [1],  mathematical  formulations  are 
developed  for  “new,  three-dimensional  packet-line  solutions  to  the  free-space 
Maxwell’s  equations.”  These  solutions  are  further  characterized  as  “real,  non¬ 
singular,  continuous  functions  which  propagate  in  a  straight  line  at  light  velocity, 


1 


They  remain  focused  for  all  time.”  It  was  subsequently  shown  by  Belanger  [2] 
that  a  whole  family  of  packet-like  solutions  of  the  homogeneous  wave  equation 
can  be  derived  that  includes  Brittingham’s  form  as  a  special  case.  Brittingham 
and  Belanger  both  observed  that  pulses  of  the  type  formulated,  while  satisfying 
Maxwell’s  equations,  carry  infinite  energy  so  that  they  are  physically  unrealizable. 
In  order  to  overcome  this  difficuly,  Brittingham  [1]  supplemented  his  three- 
dimensional  wave  packet  by  preceding  and  following  surfaces  of  discontinuity 
which  served  to  keep  the  energy  associated  with  the  packet  finite.  Work  in  this 
area  at  Harvard  began  with  the  determination  [3]  that  these  bounded  pulses  of 
Brittingham  no  longer  satisfy  Maxwell’s  equations  across  the  surfaces  of  discon¬ 
tinuity  and  are  therefore  not  physically  meaningful.  The  question  was  then  raised 
whether  there  are  other  solutions  with  the  desired  properties.  In  collaboration 
with  Professor  H.  Lehmann  of  the  University  of  Hamburg,  we  have  shown  [4] 
under  very  general  conditions  that  there  are  no  such  solutions.  In  other  words, 
the  so-called  “focus  wave  mode”  does  not  focus. 

Current  research  is  concerned  with  the  possibility  of  generating  electromag¬ 
netic  pulses  which,  while  not  “focused  for  all  time”  like  Brittingham’s,  neverthe¬ 
less  propagate  in  a  manner  such  that  the  energy  associated  with  them  decreases 
with  the  radial  distance  R  much  more  slowly  than  as  R~2.  The  justification  for 
the  belief  that  this  may  be  realizable  comes  from  recent  advances  in  particle  phy¬ 
sics.  In  particular,  it  is  learned  from  this  work  that,  even  for  extremely  short 
wavelengths,  a  logarithmic  dependence  of  such  quantities  as  the  total  cross-section 
on  the  wavelength  can,  and  does,  occur.  In  the  specific  formulation  proposed 
recently  [5],  it  is  found  that  the  energy  received  by  a  given  detector  must 
approach  zero  as  it  is  moved  further  and  further  away  from  the  pulse-generating 
antenna.  However,  this  approach  to  zero  can  be  much  slower  than  given  by  the 
form  R~2  for  an  outward-traveling  spherical  wave.  Such  cases  of  slow  decrease 


are  referred  to  as  electromagnetic  missiles.  Another  way  of  stating  the  result  is 
that  the  product  of  the  following  two  quantities  can  approach  zero  as  slowly  as 
one  wishes:  (1)  the  energy  transmitted  per  unit  area  of  receiver,  and  (2)  .he  total 
area  of  the  receiver.  Some  simple  examples  of  electromagnetic  missiles  have  been 
constructed  [5]  when  either  one,  but  not  both,  of  these  quantities  remains  finite  at 
infinite  distance.  Since  the  rate  of  decrease  with  distance  of  electromagneic  mis¬ 
siles  can  be  very  slow  and  they  have  the  advanage  of  moving  with  the  velocity  of 
light,  they  have  important  possible  applications. 

Since  the  first  example  of  an  electromagnetic  missile  in  ref.  5  involves  the 
current  distribution  on  a  circular  aperture,  a  cogent  question  is:  What  is  the  sim¬ 
plest  idealized  aperture  antenna?  From  the  theoretical  point  of  view  the  sug¬ 
gested  answer  is  the  open  end  of  a  circular  waveguide  excited  in  the  transverse 
distribution  characteristic  of  the  lowest  mode,  the  TEU  mode.  This  problem  is 
currently  being  investigated  by  Dr.  John  Myers  with  ONR  support. 
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V.  SIGNIFICANT  ACCOMPLISHMENT  REPORT 


V.l  Ultra-High-Speed  Real-Time  Optical  Processing  Via  Convolution 
and  Cross-Correlation  in  Inhomogeneously  Broadened  Optical 
Materials.  W.R.  Babbit,  Y.S.  Bai,  and  T.W.  Mossberg,  Contracts  N00014- 
84-K-0465,  ARO  DAAG-29-83-K-0040,  and  DAAG-29-84-G-0012;  Research 
Unit  6. 

The  mathematical  operation  of  convolution  (or  cross-correlation)  is  known  to 
be  important  in  various  applications  including  radar  and  pattern  recognition.  We 
have  developed  [l]  a  new  all-optical  means  of  performing  convolution  which  is  dis¬ 
tinguished  by  its  capability  of  operating  on  signals  of  extremely  large  bandwidth. 
In  this  process,  the  two  waveforms  to  be  convolved  are  temporally  encoded  on 
laser  pulses.  The  two  temporally  modulated  laser  pulses  together  with  a  third 
temporally  short  laser  pulse  pass  successively  through  an  inhomogeneously 
broadened  material  having  an  absorption  line  at  the  laser’s  frequency.  The 
material  is  thereby  stimulated  to  emit  a  delayed  optical  signal  whose  temporal 
waveform  represents  the  convolution  of  the  two  input  waveforms.  The  bandwidth 
(duration)  of  the  input  signals  is  limited  only  by  the  inhomogeneous  (inverse 
homogeneous)  absorption  bandwidth  of  the  subject  material.  In  light  of  typical 
material  parameters,  these  limits  imply  that  optical  signals  of  microsecond  dura¬ 
tion  and  multi-gigahertz  bandwidth  can  be  readily  convolved.  Acousto-optic  con¬ 
volvers  currently  in  use  [2]  cannot  match  these  bandwidths. 

In  the  last  reporting  interval,  we  have  performed  the  first  quantitative  exper¬ 
imental  study  of  this  process.  We  employed  the  convolution  process  to  perform 
mixed  binary  multiplication  at  20  MHz  bit  rates.  These  bit  rates  were  limited 
only  by  our  available  laser  modulation  equipment.  On  the  basis  of  our  measure¬ 
ments,  we  are  able  to  deduce  that  in  ideal  systems  convolution  signals  will  be 


approximately  1  percent  as  intense  as  the  input  signals  being  convolved,  and 
therefore  be  relatively  easy  to  detect. 


These  results  were  obtained  in  the  course  of  our  on-going  basic  study  of  the 
interaction  of  temporally  modulated  laser  radiation  with  inhomogeneously 
broadened  absorbers.  Such  interactions  are  important  in  a  number  of  potentially 
useful  optical  processing  and  data  storage  schemes. 
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V.2  Electromagnetic  Missiles.  T.T.Wu  and  R.W.P.  King,  Supplementary 
Research  Unit  #11A,  SOLITARY  UNATTENUATED  TRAVELING 
PULSES. 

The  SDI  has  generated  special  interest  in  the  capacity  of  electromagnetic 
pulses  to  carry  energy  from  a  high-powered  source  to  a  target.  A  newly  formu¬ 
lated  electromagnetic  wave  packet  called  a  “focus  wave  mode”  derived  by  J.N. 
Brittingham  [l] . aroused  great  interest  since  it  promised  a  permanently  focused 
three-dimensional  concentration  of  energy  that  propagates  in  a  straight  line  with 
the  speed  of  light.  The  question  was  raised:  how  can  it  be  generated? 
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A  careful  study  revealed  that  while  the  proposed  three-dimensional  wave 
packet  with  its  preceding  and  following  surfaces  of  discontinuity  satisfies 
Maxwell’s  equations,  it  does  not  satisfy  the  associated  boundary  conditions  [2].  A 
subsequent  extension  of  this  study  [3]  showed  that  no  solutions  with  the  proper- 
ti  of  the  “focus-wave  mode”  can  exist;  the  focus  waves  do  not  focus.  The 


prompt  demonstration  of  this  important  fact  is  a  significant  accomplishment. 

Does  this  mean  that  the  energy  associated  with  an  outward-traveling  elec¬ 
tromagnetic  pulse  generated  by  a  finite  antenna  must  decrease  with  distance  R 
as  l/R2—  the  limitation  on  all  continuous-wave  propagation?  A  broadly  based 
investigation  of  this  interesting  question  has  led  to  the  very  important  result  illus 
trated  by  specific  examples  that  a  decrease  in  energy  with  a  distance  that  is  much 
slower  than  l/R2  is  possible. 
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